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Abstract 

Motivated by ideas about quantum gravity, a tremendous amount of effort over the past 
decade has gone into testing Lorentz invariance in various regimes. This review summarizes 
both the theoretical frameworks for tests of Lorentz invariance and experimental advances 
that have made new high precision tests possible. The current constraints on Lorentz violat- 
ing effects from both terrestrial experiments and astrophysical observations are presented. 
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1 Introduction 



Relativity has been one of the most successful theories of the last century and is a cornerstone of 
modern physics. This review focuses on the modern experimental tests of one of the fundamental 
symmetries of relativity, Lorentz invariance. Over the last decade there has been tremendous 
interest and progress in testing Lorentz invariance. This is largely motivated by two factors. First, 
there have been theoretical suggestions that Lorentz invariance may not be an exact symmetry at 
all energies. The possibility of four dimensional Lorentz violation has been investigated in different 
quantum gravity models (including string theory [1851 111)7) . warped brane worlds |7U|. and loop 
quantum gravity 120 ]), although no quantum gravity model predicts Lorentz violation conclusively. 
Other high energy models of spacetime structure, such as non-commutative field theory, do however 
explicitly contain Lorentz violation |98j . High energy Lorentz violation can regularize field theories, 
another reason it may seem plausible. Even if broken at high energies Lorentz symmetry can still be 
an attractive infrared fixed point, thereby yielding an approximately Lorentz invariant low energy 
world [7J5] • Other ideas such as emergent gauge bosons |^ 11891 11611 180) , varying moduli [53] , 
axion-Wess-Zumino models [30], analogues of emergent gravity in condensed matter 40, 238 , ghost 
condensate [HI], space-time varying couplings [1771 15T)] . or varying speed of light cosmologies [2191 
209 also incorporate Lorentz violation. The ultimate fate of Lorentz invariance is therefore an 
important theoretical question. 

We shall primarily focus on quantum gravity induced Lorentz violation as the theoretical target 
for experimental tests. If Lorentz invariance is violated by quantum gravity, the natural scale one 
would expect it to be strongly violated at is the Planck energy of « 10 19 GeV. While perhaps 
theoretically interesting, the large energy gap between the Planck scale and the highest known 
energy particles, the trans-GZK cosmic rays of 10 11 GeV (not to mention accelerator energies of 
~ 1 TeV), precludes any direct observation of Planck scale Lorentz violation. 

Fortunately, it is very likely that strong Planck scale Lorentz violation yields a small amount of 
violation at much lower energies. If Lorentz invariance it is violated at the Planck scale, there must 
be an interpolation to the low energy, (at least nearly) Lorentz invariant world we live in. Hence 
a small amount of Lorentz violation should be present at all energies. Advances in technology 
and observational techniques have dramatically increased the precision of experimental tests, to 
the level where they can be sensitive to small low energy residual effects of Planck scale Lorentz 
violation. These experimental advances are the second factor stimulating recent interest in testing 
Lorentz invariance. One should keep in mind that low energy experiments cannot directly tell 
us whether or not quantum gravity is Lorentz invariant. Rather, they can only determine if the 
"state" that we live in is Lorentz violating. For example, it is possible that quantum gravity 
might be Lorentz invariant but contain tensor fields that acquire a vacuum expectation value at 
low energies [T85j, thereby spontaneously breaking the symmetry. Experiments carried out at low 
energies would therefore see Lorentz violation, even though it is a good symmetry of the theory at 
the Planck scale. That said, any discovery of Lorentz violation would be an important signal of 
beyond standard model physics. 

There are currently a number of different theoretical frameworks in which Lorentz symmetry 
might be modified, with a parameter space of possible modifications for each framework. Since 
many of the underlying ideas come from quantum gravity, which we know little about, the fate 
of Lorentz violation varies widely between frameworks. Most frameworks explicitly break Lorentz 
invariance, in that there is a preferred set of observers or background field other than the met- 
ric inmOS]- However others try to deform the Poincare algebra, which would lead to modified 
transformations between frames but no preferred frame (for a review see |186| L These latter 
frameworks lead to only "apparent" low energy Lorentz violation. Even further complications 
arise as some frameworks violate other symmetries, such as CPT or translation invariance, in 
conjunction with Lorentz symmetry. The fundamental status of Lorentz symmetry, broken or 
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deformed, as well as the additional symmetries makes a dramatic difference as to which experi- 
ments and observations are sensitive. Hence the primary purpose of this review is to delineate 
various frameworks for Lorentz violation and catalog which types of experiments are relevant for 
which framework. Theoretical issues relating to each framework are touched on rather briefly, but 
references to the relevant theoretical work are included. 

Tests of Lorentz invariance span atomic physics, nuclear physics, high-energy physics, relativity, 
and astrophysics. Since researchers in so many disparate fields are involved, this review is geared 
towards the non-expert /advanced graduate level, with descriptions of both theoretical frameworks 
and experimental/observational approaches. Some other useful starting points on Lorentz violation 
are |22l I27HI 11741 115.5] . The structure of this review is as follows. An general overview of various 
issues relating to the interplay of theory with experiment is given in section [5] The current 
theoretical frameworks for testing Lorentz invariance are given in sections and 0] A discussion 
of the various relevant results from earth based laboratory experiments, particle physics, and 
astrophysics is given in sections and Limits from gravitational observations are in section [7| 
Finally, the conclusions and prospects for future progress are in section |S] Throughout this review 

7/ Q/3 denotes the Minkowski (H ) metric. Greek indices will be used exclusively for spacetime 

indices whereas Roman indices will be used in various ways. Theorists' units h = c = 1 are used 
throughout. Epi denotes the (approximate) Planck energy of 10 19 GeV. 

2 General Considerations 

2.1 Defining Lorentz violation 
2.1.1 Lorentz violation in field theory 

Before we discuss Lorentz violation in general, it will be useful to detail a pedagogical example 
that will give an intuitive feel as to what "Lorentz violation" actually means. Let us work in a 
field theory framework and consider a "bimetric" action for two massless scalar fields 

S = \ j J—gd 4 x (g al3 d a <j>dpct> + (g a P + r^)d a ^) (1) 

where r Q ^ is some arbitrary symmetric tensor, not equal to g a ^ . Both g a/3 and r™' 5 are fixed back- 
ground fields. At a point, one can always choose coordinates such that g a/3 = r} a P . Now, consider 
the action of local Lorentz transformations at this point, which we define as those transformations 
for which r\ a p is invariant, on S. 1 S is a spacetime scalar, as it must be to be well-defined and 
physically meaningful. Scalars are by definition invariant under all passive diffeomorphisms (where 
one makes a coordinate transformation of every tensor in the action, background fields included) . 
A local Lorentz transformation is a subgroup of the group of general coordinate transformations so 
the action is by construction invariant under "passive" local Lorentz transformations. This implies 
that as long as our field equations are kept in tensorial form we can freely choose what frame we 
wish to calculate in. Coordinate invariance is sometimes called "observer Lorentz invariance" in 
the literature |172) although it really has nothing to do with the operational meaning of Lorentz 
symmetry as a physical symmetry of nature. 

Lorentz invariance of a physical system is based upon the idea of "active" Lorentz transforma- 
tions, where we only transform the dynamical fields <fi and ip. Consider a Lorentz transformation 
of 4> and ip, 

tf{x) = ^(A- 1 )^) (2) 

fix) = vaA- 1 )^) 

1 Note that since there are two metrics, rf^ and rf 1 ^ + T a @ there can be two different sets of transformations 
that leave one of the metrics invariant. In this sense there are two Lorentz groups. 
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where A£J is the Lorentz transformation matrix, x l{i — A^x u . The derivatives transform as 

d v <t>'{x) = (A-^d^ttA- 1 )^) (3) 
d^'(x) = (A-^d^dA- 1 )^) 

from which one can easily see that rj a ^d a (f>' (x)dp4>' (x) = rf 1 ® d a (f>(x)dp(f>(x) since by definition 
rf 1 ^ '(A _1 )^(A _1 )^ = ij^ . The ry Q/3 terms are therefore Lorentz invariant. r Q/3 is not, however, 
invariant under the action of A -1 and hence the action violates Lorentz invariance. Equations of 
motion, particle thresholds, etc. will all be different when expressed in the coordinates of relatively 
boosted or rotated observers. 

Since in order for a physical theory to be well defined the action must be a spacetime scalar, 
breaking of active Lorentz invariance is the only physically acceptable type of Lorentz violation. 
Sometimes active Lorentz invariance is referred to as "particle" Lorentz invariance |172| . We 
will only consider active Lorentz violation and so shall drop any future labelling of "observer", 
"particle" , "active" , or "passive" Lorentz invariance. For the rest of this review, Lorentz violation 
always means active Lorentz violation. For another discussion of active Lorentz symmetry in 
field theory see |24()j . Since we live in a world where Lorentz invariance is at the very least an 
excellent approximate symmetry, r"' 3 must be small in our frame. In field theoretical approaches to 
Lorentz violation, a frame in which all Lorentz violating coefficients are small is called a concordant 
frame |176| . 

2.1.2 Modified Lorentz groups 

Almost all models for Lorentz violation fall into the framework above, where there is a preferred set 
of concordant frames (although not necessarily a field theory description) . In these theories Lorentz 
invariance is broken; there is a preferred set of frames where one can experimentally determine that 
Lorentz violation is small. A significant alternative that has attracted attention is simply modifying 
the way the Lorentz group acts on physical fields. In the discussion above, it was assumed that 
everything transformed linearly under the appropriate representation of the Lorentz group. On 
top of this structure, Lorentz non-invariant tensors were introduced that manifestly broke the 
symmetry but the group action remained the same. One could instead modify the group action 
itself in some manner. A partial realization of this idea is provided by so-called "doubly special 
relativity" [151 1186) (DSR), which will be discussed more thoroughly in Sec. 13.41 In this scenario 
there is still Lorentz invariance, but the Lorentz group acts non-linearly on physical quantities. The 
new choice of group action leads to a new invariant energy scale as well as the invariant velocity 
c (hence the name doubly special). The invariant energy scale Xdsr is usually taken to be the 
Planck energy. There is no preferred class of frames in these theories, but it still leads to Lorentz 
"violating" effects. For example, there is a wavelength dependent speed of light in DSR models. 
This type of violation is really only "apparent" Lorentz violation. The reader should understand 
that it is a violation only of the usual linear Lorentz group action on physical quantities. 

2.2 Kinematics vs. dynamics 

A complete physical theory must obviously include dynamics. However, over the years a number 
of kinematic frameworks have been developed for testing Lorentz violation that do not postulate 
a complete dynamics [2461 12111 12051 12U) . Furthermore, some proposals coming from quantum 
gravity are at a stage where the low energy kinematics are partially understood/conjectured, but 
the corresponding dynamics are not understood (a good example of this is DSR 186 ). Hence 
until these models become more mature, only kinematic tests of Lorentz invariance are truly 
applicable. Strictly enforced, this rule would preclude any use of an experiment that relies on 
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particle interactions, as these interactions are determined by the dynamics of the theory. Only 
a select few observations, such as interferometry, birefringence, Doppler shifts, or time of flight 
are by construction insensitive to dynamics. However, the observational situation is often such 
that tests that use particle interactions can be applied to theories where only the kinematics is 
understood. This can be done in astrophysical threshold interactions as long as the dynamics 
are assumed to be not drastically different from Lorentz invariant physics (see section I6.4f) . In 
terrestrial experiments, one must either recognize that different experiments can give different 
values with kinematic frameworks (for an example, see the discussion of the Robertson-Mansouri- 
Scxl framework in section l3~2"|) or embed the kinematics in a fully dynamical model like the standard 
model extension fsection 14.1.1(1 . 

2.3 The role of other symmetries 

There are many other symmetries that affect how Lorentz violation might manifest itself below 
the Planck scale. The standard model in Minkowski space is invariant under four main symme- 
tries, three continuous and one discrete. There are two continuous spacetime symmetries, Lorentz 
symmetry and translation symmetry, as well as gauge and CPT symmetry. Supersymmetry can 
also have profound effects on how Lorentz violation can occur. Finally, including gravity means 
that we must take into account diffeomorphism invariance. The fate of these other symmetries in 
conjunction with Lorentz violation can often have significant observational ramifications. 

2.3.1 CPT Invariance 

Lorentz symmetry is intimately tied up with CPT symmetry in that the assumption of Lorentz 
invariance is required for the CPT theorem 162 . Lorentz violation therefore allows for (but does 
not require) CPT violation, even if the other properties of standard quantum field theory are as- 
sumed. Conversely, however, CPT violation implies Lorentz violation for local field theories |134| . 
Furthermore, many observational results are sensitive to CPT violation but not directly to Lorentz 
violation. Examples of such experiments are kaon decay fsection l5.5|) and 7-ray birefringence (sec- 
tion l6.3(l . both of which indirectly provide stringent bounds on Lorentz violation that incorporates 
CPT violation. Hence CPT tests are very important tools for constraining Lorentz violation. In 
effective field theory CPT invariance can explicitly be imposed to forbid a number of strongly 
constrained operators. For more discussion on this point see section T4. 31 

2.3.2 Supersymmetry 

SUSY, while related to Lorentz symmetry, can still be an exact symmetry even in the presence 
of Lorentz violation. Imposing exact SUSY provides another custodial symmetry that can forbid 
certain operators in Lorentz violating field theories. If, for example, exact SUSY is imposed in the 
MSSM (minimal supersymmetric standard model), then the only Lorentz violating operators that 
can appear have mass dimension five or above |137| . Of course, we do not have exact SUSY in 
nature. The size of low dimension Lorentz violating operators in a theory with Planck scale Lorentz 
violation and low energy broken SUSY has recently been analyzed in |65j . For more discussion on 
this point see Sec. 14.31 

2.3.3 Poincare Invariance 

In many astrophysics approaches to Lorentz violation, conservation of energy-momentum is used 
along with Lorentz violating dispersion relations to give rise to new particle reactions. Absence of 
these reactions then yields constraints. Energy /momentum conservation between initial and final 
particle states requires translation invariance of the underlying spacetime and the Lorentz violating 
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physics. Therefore we can apply the usual conservation laws only if the translation subgroup of the 
Poincare group is left unmodified. If Lorentz violation happens in conjunction with a modification 
of the rest of the Poincare group, then it can happen that modified conservation laws must be 
applied to threshold reactions. This is the situation in DSR: all reactions that are forbidden by 
conservation in ordinary Lorentz invariant physics are also forbidden in DSR |14fij . even though 
particle dispersion relations in DSR would naively allow new reactions. The conservation equations 
change in such a way as to compensate for the modified dispersion relations (see section l3.4l) . Due 
to this unusual (and useful) feature, DSR evades many of the constraints on effective field theory 
formulations of Lorentz violation. 

2.4 Diffeomorphism invariance and prior geometry 

If Lorentz violating effects are to be embedded in an effective field theory, then new tensors must 
be introduced that break the Lorentz symmetry (c.f. the bimetric theory of Sec. 12. If we are 
considering only special relativity, then keeping these tensors as constant is viable. However, any 
complete theory must include gravity, of course, and one should preserve as many fundamental 
principles of general relativity as possible while still introducing local Lorentz violation. There are 
three general principles in general relativity relevant to Lorentz violation: general covariance (which 
implies both passive and active diffeomorphism invariance |247) L the equivalence principle, and 
lack of prior geometry. As we saw in section|21 general covariance is automatically a property of an 
appropriately formulated Lorentz violating theory, even in flat space. The fate of the equivalence 
principle we deal with below in section ^. 51 The last principle, lack of prior geometry, is simply a 
statement that the metric is a dynamical object on the same level as any other field. Coupled with 
diffeomorphism invariance this leads to conservation of matter stress tensors (for a discussion see 
|73|L However, a fixed Lorentz violating tensor constitutes prior geometry in the same way that 
a fixed metric would. If we keep our Lorentz violating tensors as fixed objects, we immediately 
have non-conservation of stress tensors and inconsistent Einstein equations. As a specific example, 
consider again the bimetric theory We will include gravity in the usual way by adding the 
Einstcin-Hilbert Lagrangian for the metric. The resultant action is 

S = J d^x^—g (j^R + \g a0 V a 4>V^ + htf# + r^V^V^J (4) 

and the corresponding field equations are 

V a V a <j) = (5) 

V Q V Q ^ + V a (r Q/3 V^) =0 (6) 

G al3 = V Q 0V /3 - ~V a <l>V<T<t>g a f3 + - iv CT VV CT V/?a/3 - -T^VpitVv-ipgap. (7) 

Taking the divergence of J7J) and using the cf>, ip equations of motion yields 

= -V/^V^T^VerVO " \V p{T P<J V ^) (8) 

since V Q G Q ^ vanishes by virtue of the Bianchi identities. 

The right hand side of J5J) does not in general vanish for solutions to the field equations and 
therefore (JBJ) is not in general satisfied unless one restricts to very specific solutions for ip. This 
is not a useful situation, as we would like to have the full space of solutions for ip yet maintain 
energy conservation. The solution is to make all Lorentz violating tensors dynamical [1731 1155] . 
thereby removing prior geometry. If the Lorentz violating tensors are dynamical then conservation 
of the stress tensor is automatically enforced by the diffeomorphism invariance of the action. While 
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dynamical Lorentz violating tensors have a number of effects that are testable in the gravitational 
sector, most researchers have concentrated on flat space tests of Lorentz invariance where gravi- 
tational effects can be ignored. Hence for most of this review we will treat the Lorentz violating 
coefficients as fixed and neglect dynamics. The theoretical consequences of dynamical Lorentz 
violation will be analyzed only in section f4. 41 where we discuss a model of a diffeomorphism in- 
variant "aether" which has received some attention. The observational constraints on this theory 
are discussed in section 

2.5 Lorentz violation and the equivalence principle 

Lorentz violation implies a violation of the equivalence principle. Intuitively this is clear: in order 
for there to be Lorentz violation particles must travel on world-lines that are species dependent 
(and not fully determined by the mass). In various papers dealing with Lorentz violating dispersion 
relations one will sometimes see the equivalence principle being cited as a motivation for keeping 
the Lorentz violating terms equal for all particle species. We now give a pedagogical example to 
show that the equivalence principle is violated even in this case. Consider a dispersion modification 
of the form 

E 2 =m 2 +p 2 + ^\p\ 4 (9) 

for a free particle and assume is independent of particle species. If we assume Hamiltonian 
dynamics at low energy and use the energy as the Hamiltonian, then for a non-relativistic particle 
in a weak gravitational field we have 

where V(x) is the Newtonian gravitational potential m<&{x). Applying Hamilton's equations to 
solve for the acceleration yields 

to lowest order in the Lorentz violating term. From this expression it is obvious that the acceler- 
ation is mass dependent and the equivalence principle is violated (albeit slightly) for particles of 
different masses with the same f( 4 \ Of course, if the /( n) terms are different, as is natural with 
some Lorentz violating models |ll(Jj , then it is also obviously violated. As a consequence one can- 
not preserve the equivalence principle with Lorentz violation unless one also modifies Hamiltonian 
dynamics. Equivalence principle tests are therefore able to also look for Lorentz violation and vice 
versa (for an explicit example see Other examples of the relationship between equivalence 

principle violation and Lorentz violation can be found in [1401 11881 1256) . 

2.6 Systematic vs. Non-systematic violations 

Most tests of Lorentz violation deal with systematic Lorentz violation, where the deviation is con- 
stant in time/space. For example, consider the modified dispersion relation lu 2 = k 2 + f^k 4 /Ep l 
for a photon where /' 4 ' is some fixed coefficient. There is no position dependence, so the Lorentz 
violating term is a constant as the particle propagates. However, various models ',!)!)■ t2")">!l2."i II sug- 
gest that particle energy/momentum may not be constant but instead vary randomly by a small 
amount. Some authors have combined these two ideas about quantum gravity, Lorentz violation 
and stochastic fluctuations, and considered a stochastic violation of Lorentz invariance character- 
ized by a fluctuating coefficient [121 12321 1231 111)81 1115) . We will discuss non-systematic models in 
greater detail in section 1531 
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2.7 Causality and stability 
2.7.1 Causality 

It is obvious that when we introduce Lorentz violation we have to rethink causality - there is no 
universal light cone given by the metric that all fields must propagate within. Even with Lorentz 
violation we must certainly maintain some notion of causality, at least in concordant frames, since 
we know that our low energy physics is causal. Causality from a strict field theory perspective 
is usually discussed in terms of microcausality which in turn comes from the cluster decomposi- 
tion principle: physical observables at different points and equal times should be independently 
measurable. This is essentially a statement that physics is local. We now briefly review how mi- 
crocausality arises from cluster decomposition. Let 0\ (x) , O2 (y) represent two observables for a 
field theory in flat space. In a particular frame, let us choose the equal time slice t = 0, such that 
x = (0,x),y = (0, y) and further assume that x 7^ y . The cluster decomposition principle then 
states that 0\{x) and 02(2/) must be independently measurable. This in turn implies that their 
commutator must vanish, [0±(x), 02(2/)] = 0. When Lorentz invariance holds there is no preferred 
frame, so the commutator must vanish for the t — surface of any reference frame. This immedi- 
ately gives that [0\(x),02{y)\ = whenever x,y are spacelike separated, which is the statement 
of microcausality. Microcausality is related to the existence of closed timelike curves since closed 
timelike curves violate cluster decomposition for surfaces that are pierced twice by the curves. The 
existence of such a curve would lead to a breakdown of microcausality. 

Lorentz violation can induce a breakdown of microcausality, as shown in |17fi| . In this work, the 
authors find that microcausality is violated if the group velocity of any field mode is superluminal. 
Such a breakdown is to be expected, as the light cone no longer determines the causal structure and 
notions of causality based on "spacelike" separation would not be expected to hold. However, the 
breakdown of microcausality does not lead to a breakdown of cluster decomposition in a Lorentz 
violating theory, in contrast to a Lorentz invariant theory. Even if fields propagate outside the 
light cone, we can have perfectly local and causal physics in some reference frames. For example, 
in a concordant frame Lorentz violation is small, which implies that particles can be only slightly 
superluminal. In such a frame all signals are always propagated into the future, so there is no 
mechanism by which signals could be exchanged between points on the same time slice. If we 
happened to be in such a concordant frame then physics would be perfectly local and causal even 
though microcausality does not hold. 

The situation is somewhat different when we consider gravity and promote the Lorentz violating 
tensors to dynamical objects. For example in an aether theory, where Lorentz violation is described 
by a timelike four vector, the four vector can twist in such a way that local superluminal propagation 
can lead to energy- momentum flowing around closed paths |2()fi) . However, even classical general 
relativity admits solutions with closed time like curves, so it is not clear that the situation is any 
worse with Lorentz violation. Furthermore, note that in models where Lorentz violation is given 
by coupling matter fields to a non-zero, timelikc gradient of a scalar field, the scalar field also acts 
as a time function on the spacetime. In such a case, the spacetime must be stably causal (c.f. \'27'2\ ) 
and there are no closed timelike curves. This property also holds in Lorentz violating models with 
vectors if the vector in a particular solution can be written as a non- vanishing gradient of a scalar. 

Finally, wc mention that in fact many approaches to quantum gravity actually predict a failure 
of causality based on a background metric |121j as in quantum gravity the notion of a spacetime 
event is not necessarily well-defined |239) . A concrete realization of this possibility is provided in 
Bose-Einstein condensate analogs of black holes 001- Here the low energy phonon excitations obey 
Lorentz invariance and microcausality 270 . However, as one approaches a certain length scale 
(the healing length of the condensate) the background metric description breaks down and the low 
energy notion of microcausality no longer holds. 
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2.7.2 Stability 



In any realistic field theory one would like a stable ground state. With the introduction of Lorentz 
violation, one must still have some ground state. This requires that the Hamiltonian still be 
bounded from below and that perturbations around the ground state have real frequencies. It will 
again be useful to discuss stability from a field theory perspective, as this is the only framework in 
which we can speak concretely about a Hamiltonian. Consider a simple model for a massive scalar 
field in flat space similar to QJ, 



where we now assume that in some frame S the only non-zero component of r"' 3 is r 00 . This 
immediately leads to the dispersion law (1 + t 00 )E 2 = p 2 + to 2 . We can immediately deduce from 
this that if r 00 is small the energy is always positive in this frame (taking the appropriate root 
of the dispersion relation). Similar statements about energy positivity and the allowable size of 
coefficients hold in more general field theory frameworks |176) . If the energy for every mode is 
positive, then the vacuum state |0s > is stable. 

As an aside, note that while the energy is positive in S, it is not necessarily positive in a boosted 
frame S'. If t 00 > 0, then for large momentum E < p, yielding a spacelike energy momentum 
vector. This implies that the energy E' can be less than zero in a boosted frame. Specifically, for a 
given mode p in S, the energy E' of this mode in a boosted frame S' is less than zero whenever the 
relative velocity v between S and S" is greater than E/p. The main implication is that if v is large 
enough the expansion of a positive frequency mode in S in terms of the modes of S' (one can do 
this since both sets are a complete basis) may have support in the negative energy modes. The two 
vacua |0s > and |0s' > are therefore inequivalent. This is in direct analogy to the Unruh effect, 
where the Minkowski vacuum is not equivalent to the Rindler vacuum of an accelerating observer. 
With Lorentz violation even inertial observers do not necessarily agree on the vacuum. Due to 
the inequivalence of vacua an inertial detector at high velocities should see a bath of radiation just 
as an accelerated detector sees thermal Unruh radiation. A clue to what this radiation represents 
is contained in the requirement that E' < only if v > E/p, which is exactly the criteria for 
Cerenkov radiation of a mode p. In other words, the vacuum Cerenkov effect (discussed in more 
detail in section lb. 51) can be understood as an effect of inequivalent vacua. 

We now return to the question of stability. For the models in section 13.11 with higher order 
dispersion relations (E 2 — p 2 + m 2 + !^p n lE n pl 2 with n > 2) there is a stability problem for 
particles with momentum near the Planck energy if < as modes do not have positive energy 
at these high momenta. However, it is usually assumed that these modified dispersion relations are 
only effective - at the Planck scale there is a UV completion that renders the fundamental theory 
stable. Hence the instability to production of Planck energy particles is usually ignored. 

So far we have only been concerned with instability of a quantum field with a background 
Lorentz violating tensor. Dynamical Lorentz violating tensors introduce further possible instabil- 
ities. In such a dynamical theory, one needs a version of the positive energy theorem 252, 279 
that includes the Lorentz violating tensors. For aether theories, the total energy is proportional 
to the usual ADM energy of general relativity |1U4) . Unfortunately, the aether stress tensor does 
not necessarily satisfy the dominant energy condition (although it may for certain choices of co- 
efficients), so there is no proof yet that spacetimes with a dynamical aether have positive energy. 
For other models of Lorentz violation the positive energy question is completely unexplored. It is 
also possible to set limits on the coefficients of the aether theory by demanding that the theory be 
perturbatively stable, which requires that excitations of the aether field around a Lorentz violating 
vacuum expectation value have real frequencies |157| . 




(12) 
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3 Kinematic frameworks for Lorentz violation 



3.1 Systematic modified dispersion 

Perhaps the simplest kinematic framework for Lorentz violation in particle based experiments is 
to propose modified dispersion relations for particles, while keeping the usual energy-momentum 
conservation laws. This was the approach taken in much of the work using astrophysical phenomena 
in the late 1990's. In a given observer's frame in flat space, this is done by postulating that the 
usual Lorentz invariant dispersion law E 2 = p 2 + m 2 is replaced by some function E 2 = F(p,m). 
In general the preferred frame is taken to coincide with the rest frame of the cosmic microwave 
background. Since we live in an almost Lorentz invariant world (and are nearly at rest with respect 
to the CMBR), in the preferred frame F(j>,m) must reduce to the Lorentz invariant dispersion at 
small energies and momenta. Hence it is natural to expand F(p, m) about p = 0, which yields the 
expression 

E 2 = m 2 + p 2 + Fl l) p l + F\fp l p? + F i f k p l pJp k + ... (13) 

(n) 

where the constant coefficients F^ n are dimensionful and arbitrary but presumably such that the 
modification is small. The order n of the first non-zero term in depends on the underlying 
model of quantum gravity taken. Since the underlying motivation for Lorentz violation is quantum 
gravity, it is useful to factor out the Planck energy in the coefficients F^ and rewrite (|13[) as 

f (3) 

E 2 = m 2 +p 2 + E Pl f^p l + fj^pY + -S±P i P j p k + ••• (14) 

Epi 

such that the coefficients are dimensionless. 

In most of the literature a simplifying assumption is made that rotation invariance is preserved. 
In nature, we cannot have the rotation subgroup of the Lorentz group strongly broken while 
preserving boost invariance. Such a scenario leads immediately to broken rotation invariance at 
every energy which is unobserved. 2 Hence, if there is strong rotation breaking there must also be 
a broken boost subgroup. However, it is possible to have a broken boost symmetry and unbroken 
rotation symmetry. Either way, the boost subgroup must be broken. Phenomenologically, it 
therefore makes sense to look first at boost Lorentz violation and neglect any violation of rotational 
symmetry. If we make this assumption then we have 

E 2 =m 2 +p 2 + E Pl fW\p\ + f^p 2 + ^\p\ 3 + ... (15) 

tipi 

There is no a priori reason (from a phenomenological point of view) that the coefficients in 
H15|l are universal (and in fact one would expect the coefficients to be renormalized differently even 
if the fundamental Lorentz violation is universal 0) ■ We will therefore label each as 
where A represent particle species. 

Modified dispersion and effective field theory Effective field theory (EFT) is not applicable 
if one wishes to stick to straight kinematics, however the EFT implications for modified dispersion 
are so significant that they must be considered. As will be shown in detail in Sec. 14.11 universal 
dispersion relations cannot be imposed for all n from an EFT standpoint. For example, rotation- 
ally invariant n = 1,3 type dispersion cannot be imposed universally on photons |90l I23()j . The 
operators that give rise to n = 1,3 dispersion are CPT violating and induce birefringence (the 

2 In a field theory, broken rotation invariance automatically yields broken boost invariance. For example, if 
rotation invariance is broken by coupling matter to a non-zero spacelike four vector, the four vector is also not boost 
invariant. 



10 



dispersion modifications change sign based on the photon helicity). Since EFT requires different 
coefficients for particles with different properties and there is no underlying reason why all coeffi- 
cients should be the same, it is phenomenologically safest when investigating modified dispersion to 
assume that each particle has a different dispersion relation. After this general analysis is complete 
the universal case can be treated with ease. 

3.2 Robertson-Mansouri-Sexl framework 

The Robertson-Mansouri-Sexl framework [2461 12111 1277) is a well known kinematic test theory for 
parameterizing deviations from Lorentz invariance. In the RMS framework, there is assumed to 
be a preferred frame E where the speed of light is isotropic. The ordinary Lorentz transformations 
to other frames are generalized to 

t' = a- 1 (t-e-x) (16) 

3f = d- 1 x-(d- 1 -b- l )^fi-+a- l vt (17) 

v 2 

where the coefficients a, 6, d are functions of the magnitude v of the relative velocity between frames. 
This transformation is the most general one-to-one transformation that preserves rectilinear motion 
in the absence of forces. In the case of special relativity, with Einstein clock synchronization, 
these coefficients reduce to a = b^ 1 = \f\ — v 2 , d = 1. The vector e depends on the particular 
synchronization used and is arbitrary. Many experiments, such as those that measure the isotropy 
of the one way speed of light |275| or propagation of light around closed loops, have observables 
that depend on a, b, d but not on the synchronization procedure. Hence the synchronization is 
largely irrelevant and we assume Einstein synchronization. 

The RMS framework is incomplete, as it says nothing about dynamics or how given clocks and 
rods relate to fundamental particles. In particular, the coordinate transformation of Ijltjf) only has 
meaning if we identify the coordinates with the measurements made by a particular set of clocks 
and rods. If we chose a different set of clocks and rods, the transformation laws may be completely 
different. Hence it is not possible to compare the RMS parameters of two experiments that use 
physically different clocks and rods (for example, an experiment that uses a cesium atomic clock 
versus an experiment that uses a hydrogen one). However, for experiments involving a single type 
of clock/rod and light, the RMS formalism is applicable and can be used to search for violations of 
Lorentz invariance in that experiment. The RMS formalism can be made less ambiguous by placing 
it into a complete dynamical framework, such as the standard model extension of Sec. 14.1.11 In 
fact, it was shown in |179) that the RMS framework can be incorporated into the standard model 
extension. 

Most often, the RMS framework is used in situations where the velocity v is small compared 
to c. We therefore expand a, b, d in a power series in v, 

a = l + (a RMS -^)v 2 + - (18) 

b=l + (fi R MS + \)v 2 + ... (19) 
a = l + 5 RMS v 2 + ... (20) 
and will give constraints on the parameters ol R ms> PrmSt&rms instead. 

3.3 c 2 framework 

The c 2 framework |277) is the flat space limit of the THe/i |2U5| framework. This framework 
considers the motion of electromagnetically charged test particles in a spherically symmetric, static 
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gravitational field. T, H, e, \i are all parameters that fold in to the motion of the particles which 
vary depending on the underlying gravitational model. In the flat space limit, which is the c 2 
formalism, the units are chosen such that limiting speed of the test particles is one while the 
speed of light is given in terms of the THefj, parameters by c 2 = H/(Tefi). The THen and c 2 
constructions can also be expressed in terms on the standard model extension |179) . 



3.4 "Doubly special" relativity 

Doubly special relativity (DSR), which has only been extensively studied over the past few years, 
is a novel idea about the fate of Lorentz invariance. DSR is not a complete theory as it has no 
dynamics and generates problems when applied to macroscopic objects (for a discussion see :186 ). 
Furthermore, it is not fully settled yet if DSR is mathematically consistent or physically meaningful. 
Therefore it is somewhat premature to talk about robust constraints on DSR from particle threshold 
interactions or other experiments. One might then ask, why should we talk about it at all? The 
reason is twofold. First, DSR is the subject of a good amount of theoretical effort and so it is useful 
to see if it can be observationally ruled out. The second reason is purely phenomenological. As we 
shall see in the sections below, the constraints on Lorentz violation are astoundingly good in the 
effective field theory approach. With the current constraints it is difficult to fit Lorentz violation 
into an effective field theory in a manner that is theoretically natural yet observationally viable. 

DSR, if it can eventually be made mathematically consistent in its current incarnation, has 
one phenomenological advantage - it does not have a preferred frame. Therefore it evades most 
of the threshold constraints from astrophysics as well as any terrestrial experiment that looks 
for sidereal variations, while still modifying the usual action of the Lorentz group. Since these 
experiments provide almost all of the tests of Lorentz violation that we have, DSR becomes more 
phenomenologically attractive as a Lorentz violating/deforming theory. 

So what is DSR? At the level we need for phenomenology, DSR is a set of assumptions that the 
Lorentz group acts in such a way that the usual speed of light c and a new momentum scale Eosr 
are invariant. Usually Ebsr is taken to be the Planck energy - we also make this assumption. All 
we will need for this review are the Lorentz boost expressions and the conservation laws, which 
we will postulate as true in the DSR framework. For brevity we only detail the Magueijo-Smolin 
version of DSR (210) . otherwise known as DSR2 - the underlying conclusions for DSR1 ^H] remain 
the same. The DSR2 boost transformations are most easily derived from the relations 

E= i (21) 

P =T-T ( 22 ) 

where Xdsr — ^DSii> are the physical/measured energy/momentum and e,n are called the 
"pseudo-energy" and "pseudo- momentum" respectively, e, tt transform under the usual Lorentz 
transforms, which induce corresponding transformations of E and p 163 . Similarly, the e, tt for 
particles are conserved as energy and momentum normally are for a scattering problem. 3 Given 
this set of rules, for any measured particle momentum and energy, we can solve for e, tt and calculate 
interaction thresholds, etc. The invariant dispersion relation for the DSR2 boosts is given by 

E2 _ p2 = mHl-X DS nf) (23) 
1 - X DS Rin 2 



3 It may therefore seem that DSR theories are "empty", in the sense that the new definitions of E,p are merely 
mathematical manipulations without any physical meaning, i.e. e, tt represent the true energy and momentum that 
one would measure. For a discussion of this point see 12031 118rjl . Since this is not a theoretical paper, we simply 
make the DSR assumption that E,p are the measurable energy and momentum (or more specifically for analysis of 
particle interactions, the energy and momentum that are assigned to in and out states in a scattering problem). 
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This concludes our (brief) discussion of the basics of DSR. For further introductions to DSR and 
DSR phenomenology see [1861 12T1 [TH 135] , We discuss the threshold behavior of DSR theories in 
section 15.6. II 

3.5 Non-systematic dispersion 

As mentioned in Sec. 12.61 Lorentz violation is only one possibility for a signal of quantum gravity. 
Another common idea about quantum gravity is that spacetime should have a stochastic (HHl 1255) 
or "foamy" |231| structure at very small scales. Combining these two ideas has lead a number 
of authors to the idea of stochastic/non-systematic dispersion where the modifications to the dis- 
persion relation fluctuate over time. Such dispersion modifications have been phenomenologically 
parameterized by three numbers, the usual coefficient /( n ) and exponent n of Sec. 13.11 and a 
length scale L which determines the length over which the dispersion is roughly constant. After a 
particle has travelled a distance L a new coefficient /'"J is chosen based upon some model depen- 
dent probability distribution P that reflects the underlying stochasticity. Usually P is assumed 
to be a gaussian about 0, such that the average energy of the particle is given by its Lorentz 
invariant value. As well, L is in general taken to be the de Broglie wavelength of the particle in 
question. Note that in these models n is not required to be an integer, the most common choices 
are n = 5/2, n — 8/3, n = 3 |232| . 4 The only existing constraints on non-systematic dispersion 
come from threshold reactions (Sec. 16.6.2(1 and the phase coherence of light (Sec. I6.9JI . 



4 The underlying approach that yields certain choices of n has been sharply criticized in 1381 . However, the 
constraints on these models are so poor that any n is observationally feasible. 
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4 Dynamical frameworks for Lorentz violation 



4.1 Effective field theory 

The most conservative approach for a framework in which to test Lorentz violation from quantum 
gravity is that of effective field theory (EFT). Both the standard model and relativity can be 
considered EFT's and the EFT framework can easily incorporate Lorentz violation via the intro- 
duction of extra tensors. Furthermore, in many systems where the fundamental degrees of freedom 
are qualitatively different than the low energy degrees of freedom, EFT applies and give correct 
results up to some high energy scale. Hence following the usual guideline of starting with known 
physics, EFT is an obvious place to start looking for Lorentz violation. 

4.1.1 Renormalizable operators and the Standard Model Extension 

The standard model is a renormalizable field theory containing only mass dimension < 4 operators. 
If we considered the standard model plus Lorentz violating terms then we would expect a tower of 
operators with increasing mass dimension. However, without some custodial symmetry protecting 
the theory from Lorentz violating dimension < 4 operators, the lower dimension operators will 
be more important than irrelevant higher mass dimension operators (see Sec. 14.31 for details). 
Therefore the first place to look from an EFT perspective is all possible renormalizable Lorentz 
violating terms that can be added to the standard model. In |J)0 Colladay and Kostelecky derived 
just such a theory in flat space - the so-called (minimal) Standard Model Extension (mSME). 5 

One can classify the mSME terms by whether or not they are CPT odd or even. We first will 
show the terms with manifestly SU(3) x SU(2) x U(l) gauge invariance. After that, we shall give 
the coefficients in a more practical notation that exhibits broken gauge invariance. 

4.1.2 Manifestly invariant form 

We deal with CPT odd terms first. The additional Lorentz violating CPT odd operators for leptons 
are 

- {o-l) nAB~L A"f^ L b - {o-r) ^abRa!^ Rb (24) 

where La is the left-handed lepton doublet La = {va Ia)l,Ra is the right singlet (Ia)r and A, B 
are flavor indices. The coefficients (cil,r) fiAB are constant vectors that can mix flavor generations. 6 
For quarks we have similarly 

- (oq) a1 abQa7 AI( 9s - {au) nAB~U Al^U b - {o-d)iiAbDa^/^L>b (25) 
where Qa — (ha <1a)l, Ua = (ua)r, and Da = {<Ia)r- In the gauge sector we have 

(k ) K B K + {k!) K e KX ^B x B^ 
+ (k 2 ) K e KX ^Tv(W x W^ + -igWxW^W,,) 

+ (h) K e KX ^Tr(G x G^ + ^ig3G x G M G v ). (26) 

Here B^W^G^ are the U(1),SU(2),SU(3) gauge fields and B^, W^ Vl G M „ are their respective 
field strengths. The ko term in Ij26|) is usually required to vanish as it makes the theory unstable. 
The remaining a, k coefficients have mass dimension one. 

5 In the literature the mSME is often referred to as just the SME, although technically it was introduced in |90| 
as a minimal subset of an extension that involved non-renormalizable operators as well. 

6 ( a L,R)fiAB can be constant because the mSME deals with only Minkowski space. If one wishes to make the 
mSME diffcomorphism invariant, these and other coefficients would be dynamical (section 12.41 . 
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The CPT even operators for leptons in the mSME are 

^{c l )^abLaY 1 D v L b + U{c R )^ AB R Al » D v R B (27) 
while we have for quarks 

\i{c Q )»uABQ A r D v Q B + \t(cu)^A B UA^ D v U B + \i{c D )^ AB D Al » D v D B . (28) 
For gauge fields the CPT even operators are 

- \(k B ),x^B^B^ - ^(k w ) KXlll/ Tr(W KX Wn - i(fc G ) KA ^Tr(G KA G^). (29) 

The coefficients for all CPT even operators in the mSME are dimensionless. While the split 
of CPT even and odd operators in the mSME correlates with even and odd mass dimension, we 
caution the reader that this does not carry over to higher mass dimension operators. Finally, we 
will in general drop the subscripts A, B when discussing various coefficients. These terms without 
subscripts are understood to be the flavor diagonal coefficients. 

Besides the fermion and gauge field content, the mSME also has Yukawa couplings between the 
fermion fields and the Higgs. These CPT even terms are 

- \[{H l )^ ab L a ^ u Rb + {Hu)^ ab Q a ^^ v Ub + {H d )^ ab Q a ^ v Db] + h.c. (30) 
Finally, there are also additional terms for the Higgs field alone. The CPT odd term is 

iik^Y^D^ + h.c. (31) 

while the CPT even terms are 

^[{k^Y v (D^D v <t> - (k^r u 4>UB^ - (ktwytfW^] + h.c. (32) 
This concludes the description of the mSME terms with manifest gauge invariance. 



4.1.3 Practical form 

Tests of the mSME are done at low energies, where the SU(2) gauge invariance has been broken. It 
will be more useful to work in a notation where individual fermions are broken out of the doublet 
with their own Lorentz violating coefficients. With gauge breaking, the fermion Lorentz violating 
terms above give the additional CPT odd terms 

- - (33) 

and the CPT even terms 

- ^H^a^iP + iic^^y D v ^ + ^V^757 M D" i> (34) 

where the fermion spinor is denoted by if). Each possible particle species has its own set of 
coefficients. For a single particle the a M term can be absorbed by making a field redefinition 
%j) — > e~ % °"' x i\). However, in multi-particle theories involving fermion interactions one cannot re- 
move a M for all fermions |89| . However, one can always eliminate one of the a lll i.e. only the 
differences between a p for various particles are actually observable. 
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As an aside, we note that there are additional dimension < 4 U(l) invariant terms for fermions 
that could be added to the mSME once SU(2) gauge invariance is broken. These operators are 

ie„^ B'rfi- D v 1)+ '-gx^^ D v V- (35) 

These terms do not arise from gauge breaking of the renormalizable mSME in the previous sec- 
tion. However, they might arise from non-renormalizable terms in an EFT expansion. As such, 
technically they should be constrained along with everything else. However, since their origin can 
only be from higher dimension operators they are expected to be much smaller than the terms 
that come directly from the mSME. 7 

Current tests of Lorentz invariance for gauge bosons directly constrain only the electromagnetic 
sector. The Lorentz violating terms for electromagnetism are 

- \{k F ) KXtlv F KX F» v + ^{k A F) K e KXliV A x F^ (36) 

where the kp term, is CPT even and the k^F term is CPT odd. The k&F term makes the theory 
unstable, so we assume it is zero from here forward unless otherwise noted (see section lb. 31 . Now 
that we have the requisite notation to compare Lorentz violating effects directly with observation 
we turn to the most common subset of the mSME, Lorentz violating QED. 



4.1.4 Lorentz violating QED 

In many Lorentz violating tests, the relevant particles are photons and electrons, making Lorentz 
violating QED the appropriate theory. The relevant Lorentz violating operators are given by (|33l 
I34I36|) . The dispersion relation for photons will be useful when deriving birefringence constraints 
on Uf- If kp 7^ spacetime acts as a anisotropic medium, and different photon polarizations 
propagate at different speeds. The two photon polarizations, labelled e±, have the dispersion 
relation |179| 

E={l + p±a)\-f\ (37) 

where p = ^k", a 2 = ^(k a p) 2 — p 2 , k a p = (kp) aj (3sp' y P S , and p a = p a /\~pt\. Strong limits can be 
placed on this birefringent effect from astrophysical sources |179| . as detailed in Sec. 16.31 

A simplifying assumption that is often made is rotational symmetry. With rotational symmetry 
all the Lorentz violating tensors must be reducible to products of a vector field, which we denote 
by u a , that describes the preferred frame. We will normalize u a to have components (1,0,0,0) 
in the preferred frame, placing constraints on the coefficients instead. The rotationally invariant 
extra terms are 

- 6u A1 V757 M '0 + \^cu^u v i>^ D v ip + ^idUfiUu^^l^ 1 D " tjj (38) 

for electrons and 

- -^kpW^UvF^F^ (39) 

for photons. The high energy (Epi ^> E 3> m) dispersion relations for the mSME will be necessary 
later. To lowest order in the Lorentz violating coefficients they are 

E 2 = m 2 +p 2 + f^p + fPp 2 (40) 
E 2 = (1 + 4 2 V (41) 

7 Note that many of the fermion terms in the mSME can also be found in the extended Dirac equation framework 
(c.f. 1371 1193 ). Similarly, parts of the electromagnetic sector were previously known (c.f. |235l ) in the context of 
equivalence principle violations. 
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where, if s = ±1 is the helicity state of the electron, — —26s, = —(c—ds), and = fcp/2. 
The positron dispersion relation is the same as l|4U|l with the replacement p — > —p, which will change 
only the term. 

In the QED sector dimension hve operators that give rise to n — 3 type dispersion have also 
been investigated by |23(J| with the assumption of rotational symmetry. 

■J-uHF^u ■ d){u v F va ) + J-u ftffoPL + m P R )(u ■ df^ (42) 
&PI a pi 

where Pl,r are the usual left and right projection operators Pl,r = 1/2(1 ± 7 5 ) and F a @ is the 
dual of F a /3, F a P = ie a/37l5 F 7 5. One should note that these operators violate CPT. Furthermore, 
they are not the only dimension five operators, a mistake that has sometimes been made in the 
literature. For example, we could have u a u^D a D^^p. These other operators, however, do not 
give rise to n — 3 dispersion as they are CPT even. 

The birefringent dispersion relation for photons that results from (|42|) is 

f (3) E 3 

E 2 =p 2 ± tLfL. (43) 
Fpi 

for right (+) and left (-) circularly polarized photons, where f) = 2£. Similarly, the high energy 
electron dispersion is 

E 2 = m 2 +p 2 + Je( l L} (44) 
bpi 

Where ff{R,L) = 2r )R,L- 8 

We note that since the dimension five operators violate CPT, they give rise to different disper- 
sions for positrons than electrons. Whereas the coefficients for the positive and negative helicity 
states of an electron are 2ttr and 2^, the corresponding coefficients for a positron's positive 
and negative helicity states are — 2t)l and —2i]r. This will be crucially important when deriving 
constraints on these operators from photon decay. 



4.2 Non-commutative spacetime 

A common conjecture for the behavior of spacetime in quantum gravity is that the algebra of 
spacetime coordinates is actually noncommutative. This idea has led to a large amount of research 
in Lorentz violation and we would be remiss if we did not briefly discuss Lorentz violation from 
non-commutativity. We will look at only the most familiar form of spacetime non-commutativity, 
"canonical" non-commutativity, where the spacetime coordinates acquire the commutation relation 

[x a ,xp] = i—^—Q a p. (45) 

A NC 

Q a f3 is an O(l) tensor that describes the non-commutativity and A^c is the characteristic non- 
commutative energy scale. A^c is presumably near the Planck scale if the non-commutativity 
comes from quantum gravity. However, in large extra dimension scenarios Anc could be as low as 
1 TeV. For discussions of other types of non-commutativity, including those that preserve Lorentz 
invariance or lead to DSR-type theories see [1871 1225) . The phenomenology of canonical non- 
commutativity as it relates to particle physics can be found in [1471 198) . 

8 These dispersion relations also arise in some approaches to low energy dynamics from loop quantum gravity 151 
171 11201 . However, the ultimate low-energy status of Lorentz invariance in loop quantum gravity is still far from clear 
(c.f. 1505112571 1. 
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The existence of Q a p manifestly breaks Lorentz invariance and hence the size of Anc is con- 
strained by tests of Lorentz violation. However, in order to match a non-commutative theory 
to low energy observations, we must have the appropriate low energy theory, which implies that 
the infamous UV/IR mixing problem of non-commutative field theory must be tamed enough to 
create a well-defined low energy expansion. No general method for doing this is known, although 
supersymmetry |21fij can perhaps do the trick. 9 If the UV/IR mixing is present but regulated by 
a cutoff then the resulting field theory can be re-expressed in terms of the mSME [311 175) . 

In order to see how constraints come about, consider for the moment non-commutative QED. 
The Seiberg-Witten map 254 can be used to express the non-commutative fields in terms of 
ordinary gauge fields. At lowest order in A^c the effective action for low energy is then 

S = ^i^D^-m^-^F^F^ 

i e a/3 — <-> i e af3 — <-> 

A2— F *0in* D^^ + -iq T j—F a ^ 1 > 1 Dp ip 
8 4 A NC 

1 6 a P - 

4 A wc 

1 9 a/3 1 9 a/3 

-o« T3-*W. V" + ^Q-rs-FapF^F^. (46) 
A A NC 8 A NC 

Direct constraints on the dimension six non-renormalizable operators from cosmological birefrin- 
gence and atomic clocks have been considered in |7S]- A stronger bound of Ajvc > 5 TO 14 GeV |218| 
on the non-commutativity scale can be derived from clock comparison experiments with Cs/Hg 
magnetometers 0^1 ( see section 15.21 . Similarly, the possibility of constraints from synchrotron 
radiation in astrophysical systems has been analyzed in |77j . 

Other strong constraints can be derived by noting that without a custodial symmetry loop 
effects with the dimension six operators will induce lower dimension operators. In the authors 
calculated what dimension four operators would be generated, assuming that the field theory has 
some cutoff scale A. The dimension six operators induce dimension four operators of the form 
B(Q 2 ) a ^F av Fg and AQ a pQ flv F atJ -F l3v , where A, B are dimensionless numbers that depend on 
Anc, A. There are two different regimes of behavior for A, B. If A ^ Anc then A,B are 0(1) 
(up to loop factors and coupling coefficients), independent of the scale A^c- Such strong Lorentz 
violation is obviously ruled out by current experiment, implying that in this perturbative approach 
such a limit is observationally not viable. If instead one takes A << A^vc then A, B oc A 2 /A^, c . 
The resulting field theory becomes a subset of the standard model extension; specifically the new 
operators have the form of the (k F ) a [ 3l $F al3 F 1& term in It has been argued |7S| that any 

realistic non-commutative theory must eventually reduce to part of the mSME. The approach 
of |3J shows this is possible, although the presence of such a low energy cutoff must be explained. 

All of the above approaches use an expansion in @ a P , A^c to get some low energy effective 
field theory. In terms of Lorentz tests, the results are all based upon this EFT expansion and not 
on the full non-commutative theory. Therefore we will restrict ourselves to discussing limits on 
various terms in effective field theories rather than directly quoting limits on the non-commutative 
scale. We leave it up to the reader to translate this value into a constraint (if any) on Anc and or 
A. 

4.3 Symmetry and relevant /irrelevant Lorentz violating operators 

The above section illustrates a crucial issue in searches for Lorentz violation that are motivated by 
quantum gravity: why is Lorentz invariance such a good approximate symmetry at low energies? 

9 Other methods of removing UV/IR mixing exist, for an example see 12711 . 
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To illustrate the problem, let us consider the standard assumption made in much of the work on 
Lorentz violation in astrophysics - that there exist corrections to particle dispersion relations of the 
form f {n ^p n E™- 2 with n > 3 and /<") of order one. Without any protective symmetry, radiative 
corrections involving this term will generate dispersion terms of the form f( n )p 2 + Epif( n >p. These 
terms are obviously ruled out by low energy experiment. 10 Accordingly, the first place to look 
for Lorentz violation is in terrestrial experiments using the standard model extension rather than 
astrophysics with higher dimension operators. However, no evidence for such violation has been 
found. The absence of lower dimension operators implies that either there is a fine tuning in the 
Lorentz violating sector |91|. some other symmetry is present that protects the lower dimension 
operators, or Lorentz invariance is an exact symmetry. 

It is always possible that Lorentz violation is finely tuned - there are other currently unexplained 
fine-tuning problems (such as the cosmological constant) in particle physics. However, it would be 
far preferable if there was some symmetry or partial symmetry that could naturally suppress/forbid 
lower dimension operators. For rotation invariance, a discrete remnant of the original symmetry 
is enough. For example, hypercubic symmetry on a lattice is enough to forbid dimension four 
rotation breaking operators for scalars. 11 No physically meaningful equivalent construction exists 
for the full Lorentz group however (see |223| for a further discussion of this point). A discrete 
symmetry that can forbid some of the possible lower dimension operators is CPT. A number of 
the most observationally constrained operators in the mSME are CPT violating, so imposing CPT 
symmetry would explain why those operators are absent. However, the CPT even operators in the 
mSME are also very tightly bounded, so CPT cannot completely resolve the naturalness problem 
either. 

Supersymmetry is currently the only known symmetry (other than Lorentz symmetry itself) 
that can protect Lorentz violating operators of dimension four or less |137l 11591 165] . much as SUSY 
protects some lower dimension operators in non-commutative field theory (216] . If one imposes 
exact SUSY on a Lorentz violating theory, the first allowed operators are of dimension five |137) . 
These dimension five operators do not induce n = 3 type dispersion like the operators (|42|1 . Instead, 
in a rotationally invariant setting they produce dispersion relations of the form 

E 2 =p 2 +m 2 {\ + f^4- + ■■■)■ (47) 
Epi 

Such modifications are completely unobservable in astrophysical processes, although high precision 
terrestrial experiments can still probe them. Dimension 6 SUSY operators in SQED also yield 
dispersion relations that are untestable by high energy astrophysics |65] . 

Fortunately, we do not live in a SUSY world, so it may be that upon SUSY breaking appropriate 
sized operators at each mass dimension are generated. This question has recently been explored 
in For CPT violating dimension five SUSY operators in SQED, the authors find that SUSY 
breaking yields dimension three operators of the form am 2 /M, where m s is the SUSY breaking 
scale, M is the scale of Lorentz violation and a is an 0(1) coefficient. For m s as light as it could 
be (around 100 GeV), spin polarized torsion balances (see section l5~4"|) are able to place limits 
on M between 10 5 — 10 10 Ep\. It therefore is probable that these operators are observationally 
unacceptable. However, dimension five SUSY operators are CPT violating, so a combination of 
CPT invariance and SUSY would forbid Lorentz violating operators below dimension six. The 
low energy dimension four operators induced by SUSY breaking in the presence of dimension six 

10 While we are primarily concerned in this section with dimensional transmutation of higher dimension operators, 
Lorentz violating renormalizable operators for one particle of course also yield radiative corrections to other particle 
operators. For a specific example see 1291 . 

11 This can easily be seen by the following argument. Consider a kinetic term in Euclidean space for a scalar field 
<j> of the form M ab d a <f>df,<f>. In four dimensions M ab must be a dimensionless tensor that has hypercubic symmetry. 
The only such tensor is 8 ab so rotation invariance is automatically preserved. Interaction terms are by their very 
nature rotation invariant, which implies that the entire action is invariant under the full rotation group. 
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operators would then presumably be suppressed by m 2 s /M 2 . This is enough suppression to be 
compatible with current experiment if M is at the Planck scale and m s < 1 TeV. 

Another method by which Lorentz violation can occur but might have small dimension < 4 
matter operators is via extra dimension scenarios. For example, in |70| a braneworld scenario 
was considered where 4-d Lorentz invariance was preserved on the brane but broken in the bulk. 
The only particle which can then directly see Lorentz violation is the graviton - the matter fields, 
being trapped on the brane, can only feel the bulk Lorentz violation through graviton loops. The 
induced dimension < 4 operators can be quite small, depending on the exact extra-dimension 
scenario considered. Note though that this approach has been criticized in [SJ, whose authors 
argue that significant Lorentz violation in the infrared would still occur. 

In summary, the current status of Lorentz violation in EFT is mildly disconcerting for a phe- 
nomcnologist (if one really wants to believe in Lorentz violation). From an EFT point of view, 
without custodial symmetries one would expect that we would have seen signs of Lorentz violation 
by now. Imposing SUSY+CPT or a braneworld scenario may fix this problem, but then we are left 
with a model with more theoretical assumptions. Furthermore a SUSY+CPT model is unlikely to 
ever be testable with astrophysics experiments and requires significant improvement in terrestrial 
experiments to be seen |fi5) . Fortunately, since this is a phenomenological review we can blithely 
ignore the above considerations and simply classify and constrain all possible operators at each 
mass dimension. This is also the safest approach. After all, we are searching for a possible signal 
from the mysterious realm of quantum gravity and so must be careful about overly restricting our 
models. 



4.4 Lorentz violation with gravity in EFT 

The previous field theories dealt only with the possible Lorentz violating terms that can be added 
to the matter sector. Inclusion of gravity into the mix yields a number of new phenomena. Lorentz 
violating theories with a preferred frame have been studied extensively (c.f. |122l 11561 1277) and 
references therein), while an extension of the mSME into Ricmann-Cartan geometry has been 
performed in |17«ij . Ghost condensate models, in which a scalar field acquires a constant time 
derivative, thereby choosing a preferred frame, were introduced in |33|. Let us first look at the 
more generic case of |173| . 

In order to couple Lorentz violating coefficients to fermions, one must work in the vierbein 
formalism (for a discussion see |272| L In Riemann-Cartan geometry the gravitational degrees 
of freedom are the vierbein and spin connection which give the Riemann and torsion tensors in 
spacetime. For the purposes of this review we will set the torsion to zero and work strictly in 
Riemannian geometry; for the complete Lorentz violating theory with torsion see |173] (for more 
general reviews of torsion in gravity see [1431 I139*] L The low energy action involving only second 
derivatives in the metric is given by 

S = 17-^ I d4xe i R - 2A + sa0R ^ + ^R^s] (48) 

1D7TG J 

where e is the determinant of the vierbein, R, R a /3, Rap-y8 are the Ricci scalar, Ricci tensor, and 
Riemann tensor respectively, and A is the cosmological constant. G is the gravitational coupling 
constant, which can be affected by Lorentz violation. Since there is no longer translation invariance, 
in principle the Lorentz violating coefficients s"* 3 , t a ^ s vary with location, so they also behave as 
spacetime varying couplings. s a ( 3 j t af3 "' 5 can furthermore be assumed to be trace-free as the trace 
can be absorbed into G and A. There are then 19 degrees of freedom left. 

The difficulty with this formulation is that it constitutes prior geometry and generically leads 
to energy-momentum non-conservation, similar to the bimetric model in section 12.41 Again the 
matter stress tensor will not be conserved unless very restrictive conditions are placed on s a @ , t Q ^ 7<5 
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(for example that they are covariantly constant) . It is unclear whether or not such restrictions can 
be consistently imposed in a complicated metric as would describe our universe. 

A more flexible approach is to presume that the Lorentz violating coefficients are dynamical, 
as has been pursued in |122l 11561 1341 11851 1220) . In this scenario, the matter stress tensor is 
automatically conserved if all the fields are on-shell. The trade-off for this is that the coefficients 
g<*P^t a fh s must be promoted to the level of fields. In particularly they can have their own kinetic 
terms. Not surprisingly, this rapidly leads to a very complicated theory, as not only must s a/3 , t a/3lS 
have kinetic terms, but they must also have potentials that force them to be non-zero at low 
energies. (If such potentials were not present, then the vacuum state of the theory would be 
Lorentz invariant.) For generic s a P ,t a P' yS , the complete theory is not known, but a simpler theory 
of a dynamical "aether" , first looked at by |122j and expanded on by [1561 11851 111)41 151) has been 
explored. 

The aether models assume that all the Lorentz violation is provided by a vector field u a . 12 
With this assumption, s a/3 can be written as u a vfi and t a ^ s can always be reduced to an s a/3 
term due to the symmetries of the Riemann tensor. The most generic action in D-dimensions that 
is quadratic in fields is therefore 



ci, 2,3,4 are dimensionless constants, R is the Ricci scalar, and the potential V{u a u a ) is some 
function that enforces a non-zero value for u a at low energies. With a proper scaling of coefficients 
and V this value can be chosen to be unit at low energies. The model of (|49l) still allows for 
numerous possibilities. Besides the obvious choice of which coefficients are actually present, u a 
can be either space-like or time-like and in extra dimension scenarios can point in one of the four 
uncompactified dimensions or the D — 4 compactified ones. 

At low energies u a acquires an expectation value, u a , and there will be excitations Su a about 
this value. Generically, there will be a single massive excitation and three massless ones. It has 
been argued in |105j that the theory suffers stability problems unless V is of the form X(u a u a — 1) 
where A is a Lagrange multiplier. The theory is also ghost free with this potential and the further 
assumption that c\ + C4 < |133| . Assuming these conditions, aether theories possess a set 
of coupled aether-metric modes which act as new gravitational degrees of freedom that can be 
searched for with gravitational wave interferometers or by determining energy loss rates from 
inspiral systems like the binary pulsar. The same scenario generically happens for any tensor field 
that acquires a VEV dynamically (see section EJ , which implies that Lorentz violation can be 
constrained by the gravitational sector as well as by direct matter couplings. 

The aether models use a vector field to describe a preferred frame. Ghost condensate gives a 
more specific model involving a scalar field. In this scenario the scalar field <fi has a Lagrangian of 
the form P(X), where X — d a 4>d a (f>. P{X) is a polynomial in X with a minimum at some value 
X = m, i.e. (j) acquires a constant velocity at its minimum. In a cosmological setting, Hubble 
friction drives the field to this minimum, hence there is a global preferred frame determined by the 
velocity of <f>. This theory gives rise to the same Lorentz violating effects of aether theories, such 
as Cerenkov radiation and spin dependent forces |33| . In general systems that give constraints on 
the coefficients of the aether theory are likely to also yield constraints on the size of the velocity 
m. 

12 The aether models take as a starting point general relativity plus an external vector field. Recent work 11441 has 
shown that there is an alternate formulation with the same Lorentz violating consequences in terms of a non-zero 
shear field in metric-affine gravity. 




(49) 



where 
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5 Terrestrial constraints on Lorentz violation 



Having laid out the necessary theoretical background, we now discuss the various experiments and 
observations that give the best limits on Lorentz violation. 

5.1 Penning traps 

A Penning trap is a combination of static magnetic and electric fields that can keep a charged 
particle localized within the trap for extremely long periods of time (for a review of Penning traps 
see HI]). A trapped particle moves in a number of different ways. The two motions relevant for 
Lorentz violation tests are the cyclotron motion in the magnetic field and Larmor precession due 
to the spin. The ratio of the precession frequency oj s to the cyclotron frequency uj c is given by 

uj s /uj c = g/2 (51) 

where g is the g- factor of the charged particle. The energy levels for a spin 1/2 particle are given 
by = nuj c + suj s where n is an integer and s — ±1/2. For electrons and positrons, where 
j«2, the state n, s = —1/2 is almost degenerate with the state n — 1, s = +1/2. The degeneracy 
breaking is solely due to the anomalous magnetic moment of the electron and is usually denoted 
by uj a = uj s — uj c . By introducing a small oscillating magnetic field into the trap one can induce 
transitions between these almost degenerate energy states and very sensitively determine the value 
of uj a . 

The primary use of measurements of uj a is that they directly give a very accurate value of g — 2. 
However, due to their precision, these measurements also provide good tests of CPT and Lorentz 
invariance. In the mSME, the only framework that has been applied to Penning trap experiments, 
the g factor for electrons and positrons receive no corrections at lowest order. However, the 
frequencies uo a , oj c both receive corrections jSU] ■ At lowest order in the Lorentz violating coefficients 
these corrections are (with the trap's magnetic field in the z-direction) 

Vf = (1 - C g - C\ X - dyy)^° (52) 

uf = T 2b% + 2d e zo m e + 2H XY (53) 

expressed in a non-rotating frame. The unmodified frequencies are denoted by u^'S an( i the Lorentz 
violating parameters are various components of the general set given in (|3H|> and (|34(l . 

The functional form of (|52(l immediately makes clear that there are two ways to test for Lorentz 
violation. The first is to look for instantaneous CPT violation between electrons and positrons 
which occurs if the bz parameter is non-zero. The observational bound on the difference between 
u a for electrons and positrons is |w+ — uj~ \ < 2.4 ■ 10~ 21 m e [§?]■ This leads to a bound on 
bz of order bz < 10~ 21 m e . The second approach is to track uj a ^ c over time, looking for sidereal 
variations as the orientation of the experimental apparatus changes with respect to the background 
Lorentz violating tensors. This approach has been used in |217) to place a bound on the diurnal 
variation of the anomaly frequency of Aw^ < 1.6 • 10 _21 m e , which limits a particular combination 
of components of 6 M , c M „, d^, at this level. Finally , we note that similar techniques have 
been used to measure CPT violations for proton/anti-proton and hydrogen ion systems |117) . By 
measuring the cyclotron frequency over time, bounds on the cyclotron frequency variation (|52() for 
the anti-proton have established a limit at the level of 1CP 26 on components of c? . 

5.2 Clock comparison experiments 

The classic clock comparison experiments are those of Hughes [149] and Drever |1U(J| and their 
basic approach is still used today. Two "clocks", usually two atomic transition frequencies, are 
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co-located at some point in space. As the clocks move they pick out different components of the 
Lorentz violating tensors in the mSME, yielding a sidereal drift between the two clocks. The 
difference between clock frequencies can be measured over long periods, yielding extremely high 
precision limits on the amount of drift and hence the parameters in the mSME. 13 Note that this 
approach is only possible if the clocks are made of different materials or have different orientations. 

The best overall limit is in the neutron sector of the mSME and comes from a 3 He/ 129 Xe 
maser system In this setup, both noble gases are co-located. The gases are placed into 

a population inverted state by collisions with a pumped vapor of rubidium. In a magnetic field 
of 1.5 Gauss, each gas acts as a maser at frequencies of 4.9 kHz and 1.7 kHz for He and Xe, 
respectively. The Xe emission is used as a magnetometer to stabilize the magnetic field while 
the He emission frequency is tracked over time, looking for sidereal variation. At lowest order in 
Lorentz violating couplings, the Lorentz violating effect for each gas is that of a single valence 
neutron so this experiment is sensitive only to neutron parameters in the mSME. The magnitude 
of the sidereal variation A/j is given by 

27r|A/j| - | - 3.5bj + 0.012(dj - ~g DJ )\ (54) 

where J stands for the X,Y components of the Lorentz violating tensors in a non-rotating frame 
that are orthogonal to the earth's rotation axis. All parameters are understood to be the ones for 
the neutron sector of the mSME. The coefficients b, d, g are related to the mSME coefficients of 
section BH by [T75] 

bj = bj - mdjo + ^mejKLgKLO - ^jklH K l (55) 

dj = m(d 0J + d JQ ) - i(2mrf, 70 + (jklUkl) (56) 

9d,j = me,jKL(gKOL + -fi^io) - bj. (57) 

Here m is the neutron mass and cjjk is the 3-dimensional antisymmetric tensor. Barring conspir- 
atorial cancellations among the coefficients, the bound on b± = \jb\ + b Y is 6.4 ± 5.4 • 10~ 32 GeV, 
which is the strongest clock comparison limit on mSME parameters. Similarly, one can derive 
bounds on dj_,gD,± that are two to three orders of magnitude lower. Hence certain components 
of these coefficients are bounded at the level of 10 -28 GeV. A continuation of this experiment has 
recently been able to directly constrain boost violation at the level of 10~ 27 GeV J7J (sidereal 
variations look at rotation invariancc). Besides the bounds above, other clock comparison experi- 
ments |175| are able to establish the following bounds on other coefficients in the neutron sector 
of the mSME: 

\cq,j\ = \m(cj Z + c ZJ )\ < 1(T 25 GeV (58) 
|c_| = \m(c X x - c YY )\ < 1(T 27 GeV 
\cxy\ = \m(cxY + c Y x)\ < 1(T 27 GeV. 

A constraint of the dimension five operators of (|42|) for neutrons was recently derived in |52| 
using limits on the spatial variation of the hyperfine nuclear spin transition in Be + as a function of 
the angle between the spin axis and an external magnetic field |64) . Assuming the reference frame of 
the earth is not aligned with the four vector u a , the extra terms in (|42|l generically introduce a small 
orientation dependent potential into the non-relativistic Schrodinger equation for any particle. For 
Be + , the nuclear spin can be thought of as being carried by a single neutron, so this experiment 

13 Clock comparison experiments have also been used to place bounds on Lorentz violation in a conjectured low 
energy state for loop quantum gravity 12061 . 
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limits the neutron Lorcntz violating coefficients. This extra potential for the neutron leads to 
anisotropy of the hyperfine transition frequency, which can be bounded by experiment. The limits 
are roughly \r]i\ < 6 x 10~ 3 , 1 772 ] < 3 if u a is timelike and coincides with the rest frame of the 
CMBR. If u a is spacelike one has < 2 x 10~ 8 , |j7 2 | < 10~ 8 . If u a is lightlike both coefficients 
are bounded at the 10~ 8 level. Note that all these bounds are approximate, as they depend on the 
spatial orientation of the experiment with respect to spatial components of u a in the lab frame. 
The authors of [22] have assumed that the orientation is not special. 

The above constraints applies solely to the neutron sector. Other clock comparison experiments 
have been performed that yield constraints on the proton sector |84i 243 196 46 241 in the mSME. 
The best proton limit, on the b± parameter, is < 2- 10~ 27 GeV [241] . with corresponding limits 
on d±,gn.± of order 10 -25 GeV. Similar bounds have been estimated |175) from the experiment 
of Berglund et. al. |46| using the Schmidt model |251| for nuclear structure, where an individual 
nucleon is assumed to carry the entire nuclear angular momentum. The experiments of Chupp |84j . 
Prestage 243 , and Lamoreaux 196 are insensitive to proton coefficients in this model so no proton 
bounds have yet been established from these experiments. As noted in |175| . proton bounds would 
be derivable with a more detailed model of nuclear structure. 



5.3 Cavity experiments 

From the Michelson-Morley experiments onward, interferometry has been an excellent method of 
testing relativity. Modern cavity experiments extend on the ideas of interferometry and provide 
very precise tests on the bounds of certain photon parameters. The main technique of a cavity 
experiment is to detect the variation of the resonance frequency of the cavity as its orientation 
changes with respect to a stationary frequency standard. In this sense, it is similar to a clock 
comparison experiment. However, since one of the clocks involves photons, cavity experiments 
constrain the electromagnetic sector of the mSME as well. 

The analysis of cavity experiments is much easier if we make a field redefinition of the elec- 
tromagnetic sector of the mSME |179) . In analogy to the theory of dielectrics, we define two new 
fields, D,H by 

D \ = f 1 + KDE KDB \ f E 
H ) V K HE 1 + KHB J V B 

The k coefficients are related to the mSME coefficients by 



(59) 



(kde)* 1 = -2(M° j0fe (60) 

(K HB y k = ^e krs (k F )P" rs (61) 

(K DB ) jk = -{K HE y k = e k ^(k F )°^. (62) 

With this choice of fields, the modified Maxwell equations from the mSME take the suggestive 
form 

V x H - d D = (63) 

V • D = (64) 

V x E + d B = (65) 

V • B = 0. (66) 

This redefinition shows that the Lorentz violating background tensor {kp) [iVa p can be thought 
of as a dielectric medium with no charge or current density. Hence we can apply much of our 
intuition about the behavior of fields inside a dielectric to construct tests of Lorentz violation. 
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Note that since H and D depend on the components of (kF)iu>af) the properties of the dielectric 
are orientation dependent. 

Constraints from cavity experiments are not on the k parameters themselves, but rather on the 
linear combinations 



I \ll 

K tr = g( K -DBj 


(67) 


(a e +) jk = \{KDE + K HB ) jk 


(68) 


= ^(k D e - K HB ) jk - S jk k tr 


(69) 


(Ko+) jk = ^(nDB + K HE ) jk 


(70) 




(71) 



all parity even while k -\- and k - are parity odd. The usefulness of this parame- 
terization can be seen if we rewrite the Lagrangian in these parameters |179| . 

L = ^[{l + k tr )E 2 -{l-h tr )B 2 ] + ^E-[k e+ + R e -]-E-^B-^ (72) 

From this expression it is easy to see that k tr corresponds to a rotationally invariant shift in the 
speed of light. It can be shown that k e - and k + also yield a shift in the speed of light, although 
in a direction dependent manner. The coefficients K e + and k d _ control birefringence. Cavity 
experiments yield limits on « e _ and k Q+ while birefringence (section 16. 3fl bounds k e+ and k a -. 

The most straightforward way to constrain Lorentz violation with cavity resonators is to study 
the resonant frequency of a cavity. Since we have a cavity filled with an orientation dependent 
dielectric the resonant frequency will also vary with orientation. The resonant frequency of a cavity 

<*> 

where m is the mode number, c is the speed of light, n is the index of refraction (including Lorentz 
violation) of any medium in the cavity, and L is the length of the cavity. f r can be sensitive 
to Lorentz violating effects through c, n, L. Depending on the construction of the cavity some 
effects can dominate over others. For example, in sapphire cavities the change in L due to Lorentz 
violation is negligible compared to the change in c. This allows one to isolate the electromagnetic 
sector. 

In general, all cavities are sensitive to the the photon k parameters. In contrast to sapphire, 
for certain materials the strain induced on the cavity by Lorentz violation is large. This allows 
sensitivity to the electron parameters c M „ at a level equivalent to the photon parameters. Fur- 
thermore, by using a cavity with a medium, the dependence of f r on n gives additional electron 
sensitivity |226j . 

The complete bounds on the mSME coefficients for cavity experiments are given in |22U67ll2*2"71 
226 280 12071 IS2l 1261) . The strongest bounds are displayed in tabled Roughly, the components 
of K e _ and c^y are bounded at O(10~ 15 ) while k 0+ is bounded at O(10~ n ). The 10 4 difference 
arises as k 0+ enters constraints suppressed by the boost factor of the earth relative to the solar 
"rest" frame where the coefficients are taken to be constant. 

5.4 Spin polarized torsion balances 

Clock comparison experiments constrain the bj parameter for protons and neutrons. Spin polarized 
torsion balances are able to place comparable limits on the electron sector of the mSME |5Hj. The 
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Parameter 


Value (xlCT ib ) 


r e 

C XY 


0.76 ±0.35 


C YZ 


0.21 ±0.46 


C XZ 


-0.16 ±0.63 


cxx — cyy 


1.15 ±0.64 


\cxx + cyy - 2c ZZ - 0.25ref_ z | 


10 3 


Z.XY 

K e- 


-0.63 ±0.43 


K e _ 


-0.45 ±0.37 


K XZ 


0.19 ±0.37 


e x - ~< Y 


-1.3 ±0.9 




-20 ± 2 


Z.XY 

^0+ 


(0.20 ±0.21) x 10 4 


Z.YZ 
^0+ 


(0.44 ±0.46) x 10 4 


K XZ 


(-0.91 ±0.46) x 10 4 



Table 1: Cavity limits on c M „, K e _, k Q+ (Taken from |22l 12261 132l l2"61 ). Components are in a sun 
centered equatorial frame. Error bars are la. The non-zero value of Rf z is argued by the authors 
to be due to systematics in the experiment |32| . 

best limits on bi (where i is the spatial direction, including that parallel to the earth's rotation axis) 
for the electron come from two balances, one in Washington [1701 1141] and one in Taiwan |148) . 
We detail the Washington experiment for pedagogical purposes - the two approaches are similar. 
In the Washington experiment a two different types of magnets (SmCo and Alnico) are arranged in 
an octagonal shape. Four SmCo magnets are on one side of the octagon and four Alnico magnets 
are on the other. The magnetization of both types of magnets is set to be equal and in the angular 
direction around the octagon. This minimizes any magnetic interactions. However, with equal 
magnetization the net electron spin of the SmCo and Alnico magnets differs as the SmCo magnets 
have a large contribution to their overall magnetization from orbital angular momentum of Sm 
ions. Therefore the octagonal pattern of magnets has an overall spin polarization in the octagon's 
plane. 

A stack of four of these octagons are suspended from a torsion fiber in a vacuum chamber. 
The magnets give an estimated net spin polarization if equivalent to approximately 10 23 aligned 
electron spins. The whole apparatus is then mounted on a turntable. As the turntable rotates a 
bound on Lorentz violation is obtained in the following manner. Lorentz violation in the mSME 
gives rise to an interaction potential for non-relativistic electrons of the form V = bi& L , where i 
stands for direction and a 1 is the electron magnetic moment. As the turntable rotates, since b 
points in some fixed direction in space, the interaction produces a torque on the torsion balance. 
The magnet apparatus therefore twists on the torsion fiber by an amount given by 

= V H sin^o - ut)/n (74) 

where Vh is the horizontal component of V, to is the frequency of rotation, (j)o is an initial phase due 
to orientation, and k is the torsion constant. Since k and u> are known, a measurement of O will 
give the magnitude of Vh- Since a % is also known, Vh gives a limit on the size of bi. The absence 
of any extra twist limits all components of |6| for the electron to be less than 10~ 28 GeV. The 
Taiwan experiment uses a different material (DyQFe2^) |148) . The bounds from this experiment 
are of order 10 -29 GeV for the components of bi perpendicular to the spin axis and 10~ 28 GeV for 
the parallel component. 



26 



To conclude this section, we note that the torsion balance experiments are actually sensitive 
enough to also constrain the dimension 5 operators in l|42[l . Assuming that all lower dimension 
operators are absent, the constraint on the dimension five operators is \t)r — r\i\ < 4 230 . 



5.5 Neutral mesons 

Mesons have long been used to probe CPT violation in the standard model. In the framework of 
the mSME, CPT violation also implies Lorentz violation. Let us focus on kaon tests, where most 
of the work has been done. The approach for the other mesons is similar 169; |J. The relevant 
parameter for CPT and Lorentz violation in neutral kaon systems is for the down and strange 
quarks (since K = ds). As we mentioned previously, one of the can always be absorbed by a 
field redefinition. Therefore only the difference between the quark a M 's, Aa M = r^cc — r s afj controls 
the amount of CPT violation and is physically measurable. Here r^s are coefficients that allow for 
effects due to the quark bound state 



A generic kaon state ^>k is a linear combination of the strong eigenstates, K°, K . If we write 
^ k in two component form the time evolution of the ^> k wavefunction is given by a Schrodinger 
equation: 

Kf°°)=K?) (75) 

where the Hamiltonian H is a 2 x 2 complex matrix. H can be decomposed into real and imaginary 
parts, H = M — iF. M and T are hermitian matrices usually called the mass matrix and decay 
matrix respectively. The eigenstates of H are the physically propagating states, which are the 
familiar short and long decay states K$ and K^. CPT violation only occurs when the diagonal 
components of H are not equal |198| . In the mSME, the lowest order contribution to the diagonal 
components of H occurs in the mass matrix M, contributions to T are higher order |184| . Hence 

the relevant observable for this type of CPT violation in the kaon system is the K°,K mass 
difference, Ak = \mK a ~ m -^°\/ rn K°- 14 

In the mSME the deviation Ak is (as usual) orientation dependent. In terms of Aa M , we 
have EU 

A*«|^| (76) 
m K o 

where /3 M is the four- velocity of the kaon in the observer's frame. The mass difference Ak has been 
extremely well measured by experiments such as KTeV |234) or FNAL E773 at Fermilab 253]. By 
looking for sidereal variations or other orientation effects one can derive bounds on each component 
of a^. The best current bounds to not quite achieve this but rather constrain a combination of 
parameters. A linear combination of Aao, Aa z is bounded at the level of 10~ 20 GeV |168| and a 
combination of Aa x , Aa v is constrained at the 10~ 21 GeV 12341 level. 



5.6 Doppler shift of lithium 

If Lorentz invariance is violated then the transformation laws for clocks with relative velocity will 
be different from the usual time dilation. The RMS framework of section 13.21 provides a conve- 
nient parameterization of how the Doppler shift can deviate from its standard relativistic form. 
Comparisons of oscillator frequencies under boosts therefore can constrain the ctRMS parameter in 
the RMS framework. The best test to date comes from spectroscopy of lithium ions in a storage 
ring 249 . In this experiment, "' Li + ions are trapped in a storage ring at a velocity of 0.064c. 
The transition frequencies of the boosted ions are then measured and compared to the transition 
frequencies at rest, providing a bound on the deviation from the special relativistic Doppler shift 

14 For a more thorough discussion of CPT (and CP) tests, see for example 11991 . 
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of |orms| < 2 • 10 7 in the RMS framework. Recently, the results of |249] have been reinterpreted 
in the context of the mSME. For the electron/proton sector the approximate bounds are (197) 

\c? xx + <? YY -2<? zz \<lQ- 11 

\<? tj + ( ? jt \<iq-» 

\<?xx+<5rY-Mzz\ <!0" 5 

\c e TJ + <JT\ <1(T 2 (77) 

where J = X, Y, Z in a heliocentric frame. In the photon sector, the limit ktr < O(10 -5 ) can also 
be set from this experiment |2fi9| . 

5.7 Muon experiments 

Muon experiments provide another window into the lepton sector of the mSME. As discussed in 
section EDfl if the mSME coefficients are to be small then there must be some small energy scale 
suppressing the Lorentz violating coefficients. There are only a few available small scales, namely 
particle masses or a symmetry breaking scale. If we assume the scale is particle mass, then muon 
based experiments would have a signal at least 10 2 larger than equivalent electron experiments due 
to the larger mass of the muon. The trade-off, of course, is that muons are unstable so experiments 
are intrinsically more difficult. 

There are two primary experiments that give constraints on the muon sector. First, spin 
transitions in muonium (/i + e _ ) have been used to place a bound on bj for the muon (see (|55|) for 
the definition of bj) |150| . Even though muonium is a muon-electron system, the muon sector of 
the mSME can be isolated by placing the muonium in a strong magnetic field and looking for a 
particular frequency resonance that corresponds to muon spin flips. The sidereal variation of this 
transition frequency is then tracked yielding a limit on bj of 

\bj\ < 5 x KT 22 77v (78) 

where J = X, Y in a non-rotating frame with Z oriented along the earth's spin axis. 

The second muon experiment that yields strong limits is the \i~ /fj, + g-2 experiment [5511551172"] . 
In this experiment relativistic \i~ (or n + ) are injected into a storage ring and allowed to decay. 
The deposit rate of the decay products along the detector is sensitive to the evolution of the spin 
of the muon, which in turn is a function of g — 2 for the muon. Lorentz violation changes this 
evolution equation, and therefore this type of g — 2 experiments can bound the mSME. As in the 
case of the g — 2 experiments in section l5.ll two types of bounds can be placed from the muon g — 2 
experiment. The first is a direct comparison between the g— 2 factors for jjT and which limits 
the CPT violating coefficient bz < 10~ 22 GeV. Furthermore, an analysis of sidereal variations 
involving only one of the [iT / /i + at the current sensitivity in |72) could bound the b j coefficient 
at the level of 1CT 25 GeV [HE]- 

5.8 Constraints on the Higgs sector 

Since the constraints on various parameters of the mSME are so tight, one can derive interesting 
indirect constraints on unmeasured sectors by considering loop effects. Such an approach has been 
recently taken in [28], where loop corrections to mSME coefficients from Lorentz violation in the 
Higgs sector are considered. Such an approach could be used with any particle, but since the Higgs 
is an observationally hidden sector such an analysis is more important as direct tests are unlikely 
any time soon. There are four parameters in the Higgs sector of the mSME (see section l4.1.1|l . 

Constraints on the antisymmetric part of (k^Y" , which we denote (k^)^", and (k^sY") (k^wY' 
come from the birefringence constraints on photon propagation (see section lf)3|) . Here the loop cor- 
rections to the photon propagator induce a non-zero (kp)^^, which can be directly constrained. 
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This yields a constraint on all three coefficients of order 10~ 16 . A bound {k^Y" < 10~ 13 can be 
derived from the cyclotron frequencies of hydrogen ions and anti-protons. Bounds on the CPT 
violating term {k,pY come from both the spin polarized torsion balance experiments and the noble 
gas maser. The torsion experiments bound the t, z components (where z is parallel to the earth's 
rotation axis) at the level of 10~ 27 GeV and the transverse components at 10~ 25 GeV. The He/Xe 
maser system gives a better, although less clean, bound on the transverse components of order 
10" 31 GeV. 
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6 Astrophysical constraints on Lorentz violation 



6.1 Relevance of astrophysical observations 

Terrestrial experiments are invariably concerned with low energy processes. They are therefore 
best suited for looking at the mSME, which involves lower dimension operators. Astrophysics is 
more suited for directly constraining higher dimension operators as the Lorentz violating effects 
scale with energy. As mentioned in section 14.31 the existence of Lorentz violating higher dimen- 
sional operators would generically generate lower dimension ones. At the level of sensitivity of 
astrophysical tests, the size of the corresponding lower dimension operators should give signals in 
terrestrial experiments. Hence if a signal is seen in astrophysics for Lorentz violation, one must 
then explain why Lorentz invariance passes all the low energy tests. As mentioned in section T4. 31 
exact SUSY, which is the only known mechanism to completely protect lower dimension operators, 
yields dispersion modifications (the primary method used in astrophysics) that are unobservable. 
In summary there is currently no "natural" and complete way that astrophysics might observe 
Lorentz violation but terrestrial experiments confirm Lorentz invariance. That said, physics is 
often surprising, and it is therefore still important to check for Lorentz violating signals in all 
possible observational areas. 



6.2 Time of flight 

The simplest astrophysical observations that provide interesting constraints on Planck scale Lorentz 
violation are time of flight measurements of photons from distant sources |2()l HOfil 1109) . This is 
also one of two processes (the other being birefringence) that can be directly applied to kinematic 
models. With a modified dispersion relation of the form l|15fl and the assumption that the velocity 
is given by v = dE/dp 15 , the velocity of a photon is given by 

If n ^ 2 the velocity is a function of energy, and the time of arrival difference AT between two 
photons at different energies travelling over a time T is 

AT = AvT = - Ul v 1 „ 2 - — -T 80 

2E n p f 

where E\^ are the photon energies. The large time T plays the role of an amplifier in this process, 
compensating for the small ratio E/Epi. 16 For n = 1 there are much better low energy constraints, 
while for n = 4 the constraints are far too weak to be useful. Hence we shall concentrate on n = 3 
type dispersion, where this constraint has been most often applied. 

The best limits |53j are provided by observations of rapid flares from Markarian 421, a blazar 
at a redshift of approximately z = 0.03, although a number of other objects give comparable 
results |250l 162) . The most rapid flare from Markarian 421 showed a strong correlation of flux at 

15 In K-Minkowski space there is currently some debate as to whether the standard relation for group velocity is 
correct 19 11861 . Until this is resolved, v = dE/dp remains an assumption that might be modified in a DSR context. 
It obviously holds in field theoretic approaches to Lorentz violation. 

le This is the first example of a significant constraint on terms in particle dispersion/effective field theory that are 
Planck suppressed, which would naively seem impossible. The key feature of this reaction is the interplay between 
the long travel time and the large Planck energy. In general any experiment that is sensitive to Planck suppressed 
operators is either extremely precise (as in terrestrial tests of the mSME) or has some sort of "amplifier". An 
amplifier is some other scale (such as travel time or particle mass) which combines with the Planck scale to magnify 
the effect. 
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1 and 2 TeV on a timescale of 280 seconds. If we assume that the flare was emitted from the same 
event at the source, the time of arrival delay between 1 and 2 TeV photons must be less than 280 
s. Combining all these factors yields the limit \f^\ < 128. 

A possible problem with the above bound is that in a single emission event it is not known 
if the photons of different energies are produced simultaneously. If different energies are emitted 
at different times that might mask a LV signal. One way around this is to look for correlations 
between time delay and redshift, which has been done for a set of GRB's in |1()9) . Since time of 
flight delay is a propagation effect that increases over time, a survey of GRB's at different redshifts 
can separate this from intrinsic source effects. This enables constraints to be imposed (or LV to be 
observed) despite uncertainty regarding source effects. The current data from GRB's limit /( 3 ) to 
be less than O(10 3 ) |109| . Therefore significant observational progress must be made in order to 
reach 0(1) bounds on Z^ 3 -*. Improvements on this limit might come from observations of gamma 
ray bursts with new instruments such as GLAST, however concerns have been raised that source 
effects may severely impair this approach as well |242l 1106) . Higher order dispersion corrections 
seem unlikely to ever be probed with time of flight measurements. 

The limit |/^ 3) | < 128 can be easily applied to the EFT operators in lj32jl. From (gSJ we 
trivially see that the constraint on £ is |£| < 64, again comparing the 2 TeV peak to the 1 TeV. 
It might seem that we can get a better constraint by demanding the time delay between 2 TeV 
right handed and left handed photons is less than 280 seconds. However, the polarization of the 
flare is unknown, so it is possible (although perhaps unlikely) that only one polarization is being 
produced. If one can show that both polarizations are present, then one can further improve this 
constraint. However, the time of flight constraints are much weaker than other constraints that 
can be derived on the operators in l)42|l from birefringence so this line of research would not be 
fruitful. 

DSR theories may also predict a time of flight signal, where the speed of light is effectively 
given by the group velocity of an n — 3 type dispersion relation. 17 If there is such a frequency 
dependence, it is not expected that DSR also yields birefringence as in the EFT case. An n = 3 
type dispersion for photons without birefringence would hence be a strong signal for DSR or 
something similar. Coupled with the fact that DSR does not affect threshold reactions or exhibit 
sidereal effects, time of flight analyses provide the only currently realistic probe of DSR theories. 
Unfortunately, since the invariant energy scale is usually taken to be the Planck energy, time of 
flight constraints are still one to two orders of magnitude below what is needed to constrain/probe 
DSR. 

As an aside, note that the actual measurement of the dependence of the speed of light with fre- 
quency in a telescope such as GLAST L26D1 has a few subtleties in a DSR framework. Let us make 
the (unrealistic) assumption that the situation is as good as it could possibly be experimentally: 
there is a short, high energy gamma ray burst from some astrophysical source where all the photons 
are emitted from the same point at the same time. The expected observational signal is then a 
correlation between the photon time of arrival and energy. The time of arrival is fairly straight- 
forward to measure, but the reconstruction of the initial photon energy is not so easy. GLAST 
measures the initial photon energy by calorimctry - the photon goes through a conversion foil and 
converts to an electron-positron pair. The pair then enters a calorimeter, which measures the en- 
ergy by scintillation. The initial particle energy is then only known by reconstruction from many 
events. Energy reconstruction requires addition of the multitude of low energy signals back into 
the single high energy incoming photon. Usually this addition in energy is linear (with corrections 
due to systematics/experimental error). However, if we take the DSR energy summation rules as 
currently postulated the energies of the low energy events add non-linearly, leading to a modified 
high energy signal. One might guess that since the initial particle energy is well below the Planck 

17 There currently seems to be some disagreement about this. For example, in 21 258 an energy dependent speed 
of light is argued for, whereas in 186 no such modification is found. 
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scale the non-linear corrections make little difference to the energy reconstruction. However, to 
concretely answer such a question the multi-particle sector of DSR must be properly understood 
(for a discussion of the problems with multi-particle states in DSR see |186p . 

Finally, while photons are the most commonly used particle in time of flight tests other particles 
may also be employed. For example, it has been proposed in [HJ that neutrino emission from GRB's 
may also be used to set limits on n = 3 dispersion. Observed neutrino energies can be much higher 
than the TeV scale used for photon measurements, hence one expects that any time delay is greatly 
magnified. Neutrino time delay might therefore be a very precise probe of even n > 3 dispersion 
corrections. Of course, first an identifiable GRB neutrino flux must be detected, which has not 
happened yet jjp. Assuming that a flux is seen and able to be correlated on the sky with a GRB, 
one must still disentangle the signal. In a DSR scenario, where time delay scales uniformly with 
energy this is not problematic, at least theoretically. However, in an EFT scenario there can be 
independent coefficients for each helicity, thereby possibly masking an energy dependent signal. 
For n — 3 this complication is irrelevant if one assumes that all the neutrinos are left-handed (as 
would be expected if produced from a standard model interaction) as only would then apply. 
For n > 3 the possible operators are not yet known, so it is not clear what bounds would be set 
by limits on neutrino time of flight delays. 

6.3 Birefringence 

A constraint related to time of flight is birefringence. The dimension five operators in (|42[1 as 
well as certain operators in the mSME induce birefringence - different speeds for different photon 
polarizations (|43() 18 A number of distant astrophysical objects exhibit strong linear polarization in 
various low energy bands (see for example the sources in [1781 1127] ). Recently linear polarization 
at high energies from a GRB has been reported [85] . though this claim has been challenged |248l 
1274) . Lorentz violating birefringence can erase linear polarization as a wave propagates, hence 
measurements of polarization constrain the relevant operators. 

The logic is as follows. We assume for simplicity the framework of Sec. 14.1.41 and rotation 
invariance; the corresponding analysis for the general mSME case can be found in |178| . At the 
source, assume the emitted radiation is completely linearly polarized, which will provide the most 
conservative constraint. To evolve the wave, we must first decompose the linear polarization into 
the propagating circularly polarized states, 

A a = e a L e~ ikL - x + e a R e~ ikR - x . (81) 

If we choose coordinates such that the wave is travelling in the z-direction with initial polarization 
(0,1,0,0) then €l = (0,1, — i, 0) and en = (0, 1, i, 0). Rearranging slightly we have 

A" = e ~^Lt-k z z) ( e a + £ « e -«( UJ -^)() ( 82 ) 

which describes a wave with a rotating polarization vector. Hence in the presence of birefringence 
a linearly polarized wave rotates its direction of polarization during propagation. This fact alone 
has been used to constrain the Uaf term in the mSME to the level of 10~ 42 GeV by analyzing the 
plane of polarization of distant galaxies |74j . 

A variation on this constraint can be derived by considering birefringence when the difference 
Au) = ujji — lul is a function of k z . A realistic polarization measurement is an aggregate of the 
polarization of received photons in a narrow energy band. If there is significant power across the 
entire band, then a polarized signal must have the polarization direction at the top nearly the 

18 Gravitational birefringence has also been studied extensively in the context of non-metric theories of gravitation, 
which also exhibit Lorentz violation. See for example 259 116 244 for discussions of these theories and the parallels 
with the mSME. 
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same as the direction at the bottom. If the birefringence effect is energy dependent, however, 
the polarization vectors across the band rotate differently with energy. This causes polarization 
"diffusion" as the photons propagate. Given enough time the spread in angle of the polarization 
vectors becomes comparable to 2ir and the initial linear polarization is lost. Measurement of linear 
polarization from distant sources therefore constrains the size of this effect and hence the Lorentz 
violating coefficients. We can easily estimate the constraint from this effect by looking at when 
the polarization at two different energies (representing the top and bottom of some experimental 
band) is orthogonal, i.e A" Et ^A^ Eb ^ = 0. Using (|52|) for the polarization gives 

\&U{E T ) ~ Aw(£ B )|i ps 7r. (83) 

Three main results have been derived using this approach. Birefringence has been applied to 
the mSME in [1781 1179) . Here, the ten independent components of the two coefficients k e+ and 
k _ (see section that control birefringence are expressed in terms of a ten-dimensional vector 
k a |179| . The actual bound, calculated from the observed polarization of sixteen astrophysical 
objects, is \k a \ < 10~ 32 . 19 A similar energy band was used to constrain £ in (|42() to be |£| < 
O(10 -4 ) |127j . Recently, the reported polarization of GRB021206 jHS] was used to constrain £ to 
|£| < O(10~ 14 ) |158| . but since the polarization claim is uncertain [2481 1274] such a figure cannot 
be treated as an actual constraint. 



6.4 Threshold constraints 

We now turn our attention from astrophysical tests involving a single particle species to threshold 
reactions, which often involve many particle types. Before delving into the calculational details of 
energy thresholds with Lorentz violation, we give a pedagogical example that shows why particle 
decay processes (which involve rates) give constraints that are only functions of reaction threshold 
energies. Consider photon decay, 7 — > e + e~ (see section fo.5.51 for details). In ordinary Lorentz 
invariant physics the photon is stable to this decay process. What forbids this reaction is solely 
energy/momentum conservation - two timelike four-momenta (the outgoing pair) cannot add up 
to the null four momentum of the photon. If, however, we break Lorentz invariance and assume a 
photon obeys a dispersion relation of the form 

- 2 = fc 2 + /f^ (84) 

while electrons/positrons have their usual Lorentz invariant dispersion, then it is possible to satisfy 
energy conservation equation if > (to see this intuitively, note that the extra term at high 
energies acts as a large effective mass for a photon). Therefore a photon can decay to an electron 
positron pair. 

This type of reaction is called a threshold reaction as it can happen only above some threshold 
energy uj t h ~ (rn^Epi/ f^) 1 / 3 where m e is the electron mass. The threshold energy is translated 

(3) 

into a constraint on fj in the following manner. We see 50 TeV photons from the Crab neb- 



ula |268) . hence this reaction must not occur for photons up to this energy as they travel to us 
from the Crab. If the decay rate is high enough, one could demand that u t h is above 50 TeV, 
constraining an d limiting this type of Lorentz violation. For 0(1) , u>th ~ 10 TeV, and so 
we can get a slightly better than 0(1) constraint on /^ 3 ' from 50 TeV photons 152 . If, however, 
the rate is very small then even though a photon is above threshold it could still reach us from 

19 This severe limit on birefringence provides an interesting limitation to allowable spacetime metrics in the ap- 
proach of Hehl and others 142 . In this approach, linear constitutive relations for electromagnetism are postulated 
as fundamental and the metric is derived from the constitutive relation. A lack of birefringence implies that the 
metric must be Riemannian in this approach 11941 . 
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the Crab. Using the Lorentz invariant expression for the matrix element M (i.e. just looking at 
the kinematical aspect of Lorentz violation) one finds that as ui increases above u>th the rate very 
rapidly becomes proportional to lo 2 / Epi. If a 50 TeV photon is above threshold the decay time 
is then approximately 10 _11 //.y seconds. The travel time of a photon from the Crab is ~ 10 11 
seconds. Hence if a photon is at all above threshold it will decay almost instantly relative to the 
observationally required lifetime. Therefore we can neglect the actual rate and derive constraints 
simply by requiring that the threshold itself is above 50 TeV. 

It has been argued that technically, threshold constraints can't truly be applicable to a kine- 
matic model where just modified dispersion is postulated and the dynamics/matrix elements are 
not known. This isn't actually a concern for most threshold constraints. For example, if we wish to 
constrain at 0(1) by photon decay then we can do so as long as M. is within 11 orders of mag- 
nitude of its Lorentz invariant value (since the decay rate goes as |A4| 2 ). Hence for rapid reactions, 
even an enormous change in the dynamics is irrelevant for deriving a kinematic constraint. Since 
kinematic estimates of reaction rates are usually fairly accurate (for an example see [2021 [201 ) 
one can derive constraints using only kinematic models. In general, under the assumption that 
the dynamics is not drastically different from that of Lorentz invariant effective field theory, one 
can effectively apply particle reaction constraints to kinematic theories since the decay times are 
extremely short above threshold. 

There are a few exceptions where the rate is important, as the decay time is closer to the travel 
time of the observed particle. Any type of reaction involving a weakly interacting particle such 
as a neutrino or graviton will be far more sensitive to changes in the rate. For these particles, 
the decay time of observed particles can be comparable to their travel time. As well, any process 
involving scattering, such as the GZK reaction (p + jcmbr — ► P + or photon annihilation 
(27 — > e + + e~) is more susceptible to changes in M. as the interaction time is again closer to the 
particle travel time. Even for scattering reactions, however, M. would need to change significantly 
to have any effect. Finally, M is important in reactions like (7 — > 37), which are not observed in 
nature but do not have thresholds [1541 11831 l3l l2"l IT^l] . In these situations, the small reaction rate 
is what may prevent the reaction from happening on the relevant timescales. For all of these cases, 
kinematics only models should be applied with extreme care. We now turn to the calculation of 
threshold constraints assuming EFT. 

6.5 Particle threshold interactions in EFT 

When Lorentz invariance is broken there are a number of changes that can occur with threshold 
reactions. These changes include shifting existing reaction thresholds in energy, adding additional 
thresholds to existing reactions, introducing new reactions entirely, and changing the kinematic 
configuration at threshold [86111301115411200] . By demanding that the energy of these thresholds is 
inside or outside a certain range (so as to be compatible with observation) one can derive stringent 
constraints on Lorentz violation. 

In this section we will describe various threshold phenomena introduced by Lorentz violation 
in EFT and the constraints that result from high energy astrophysics. Thresholds in other models 
are discussed in section We will use rotationally invariant QED as the prime example when 
analyzing new threshold behavior. The same methodology can easily be transferred to other 
particles and interactions. A diagram of the necessary elements for threshold constraints and the 
appropriate sections of this review is shown in Figure ^ 

Thresholds are determined by energy-momentum conservation. Since we are working in straight 
EFT in Minkowski space, translational invariance implies that the usual conservation laws hold, 
i.e. + + ■■■ = Pa + Pa + •••) where p a is the four momentum of the various particles 
A, B,C, D.... Since this just involves particle dispersion, we can neglect the underlying EFT for 
the general derivations of thresholds and threshold theorems. EFT comes back into the picture 
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Assumptions (6.5.2) 

a) Rotation Invariance 

b) Dispersion Monotonicity 

c) Energetic incoming particle 



Threshold Theorems (6.5.3) 

a) Outgoing particles parallel 

b) Initial particles anti-parallel 

c) Minimum energy final state 



Dispersion (3.1) 

E l =p* + m 1 + f™ 
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New Threshold Phenomena (6.5.4) 

a) Asymmetric thresholds 

b) Hard Cerenkov thresholds 

c) Upper thresholds 
dIHelicity decay 



Constraints in QED (6.5.5) 

a) Photon decay 

b) Vacuum Cerenkov 

| c)Photon annihilation 

Constraints from UHECR (6.5.6) 
a)GZK cutoff 

bQUHECR Cerenkov 

Constraints from Gravitational 
Cerenkov (6.5.7) 



Figure 1: Elements involved in threshold constraints. 



when we need to determine a) the actual dispersion relations that occur in a physical system to 
establish constraints and b) matrix elements for actual reaction rates (c.f. |2U1| ). 

Threshold constraints have been looked at for reactions which have the same interaction vertices 
as in Lorentz invariant physics. The reaction rate is therefore suppressed only by gauge couplings 
and phase space, n > 2 dispersion requires higher mass dimension operators, and these opera- 
tors will generically give rise to new interactions when the derivatives are made gauge covariant. 
However, the effective coupling for such interactions is the same size as the Lorentz violation and 
hence is presumably very small. These reactions are therefore suppressed relative to the Lorentz 
invariant coupling and can most likely be ignored, although no detailed study has been done. 



6.5.1 Required particle energy for "Planck scale" constraints 

We now give another simple example of constraints from a threshold reaction to illustrate the 
required energy scales for constraints on Planck scale Lorentz violation. The key concept for 
understanding how threshold reactions are useful is that, as we briefly saw for the photon decay 
reaction in section ing! particle thresholds arc determined by particle mass, which is a small number 
that can offset the large Planck energy. To see this in more detail, let us consider the vacuum 
Cerenkov effect, A —> A + 7, where A is some massive charged particle. In usual Lorentz invariant 
physics, this reaction does not happen due to energy-momentum conservation. However, consider 
now a Lorentz violating dispersion relation for A of the form 

E 2 =p 2 + m 2 + f^^ (85) 
^pi 

(n) 

with f A > 0. For simplicity in this pedagogical example we shall not change the photon dispersion 
relation u) — k. Cerenkov radiation usually occurs when the speed of the source particle exceeds 
the speed of light in a medium. The same analysis can be applied in this case, although for more 
general Lorentz violation there are other scenarios where Cerenkov radiation occurs even though 
the speed condition is not met (see below) |154j . The group velocity of A, v = dE/dp, is equal to 
one at a momentum 



Pth 



(n-l)/i n) 



l/n 

(86) 



I A 

and so we see that the threshold momenta can actually be far below the Planck energy, as it is 
controlled by the particle mass as well. For example, electrons would be unstable with n = 3 
and f y e > of 0(1) at 10 TeV, well below the maximum electron energies in astrophysical systems. 
We can rewrite <|86[1 for and see that the expected constraint from the stability of A at some 
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momentum p b s is 

(» " l )P obs 

Therefore constraints can be much less than order one with particle energies much less than Epi. 
The orders of magnitude of constraints on estimated from the threshold equation alone (i.e. 
we have neglected the possibility that the matrix elements are small) for various particles are given 
in Table H 152. 



Table 2: Orders of magnitude of vacuum Cerenkov constraint for various particles 



Particle 


V 


e 


P + 


mass 


< 1 eV 


0.511 MeV 


938 GeV 


Fobs 


100 TcV 


50 TeV 


10 20 eV 


n = 2 


1Q -28 


10~ 15 


1Q -22 


n = 3 


1Q -14 


10- 2 


10- 14 


n — 4 


1 


10 12 


10~ 6 



For neutrinos, p a i, s comes from AMANDA data |123) . The p b s for electrons comes from the 
expected energy of the electrons responsible for the creation of ~ 50 TeV gamma rays via inverse 
Compton scattering 188 268 in the Crab nebula. For nrotons. the is from AGASA data |267) . 

We include the neutrino, even though it is neutral, since neutrinos still have a non-vanishing 
interaction amplitude with photons. We shall talk more about neutrinos in section 16.81 The 
neutrino energies in this table are those currently observed; if future neutrino observatories see 
PeV neutrinos (as expected) then the constraints will increase dramatically. 

This example is overly simplified, as we have ignored Lorentz violation for the photon. However, 
the main point remains valid with more complicated forms of Lorentz violation: constraints can 
be derived with current data that are much less than 0(1) even for n — 4 Lorentz violation. We 
now turn to a discussion of the necessary steps for deriving threshold constraints, as well as the 
constraints themselves for more general models. 

6.5.2 Assumptions 

One must make a number of assumptions before one can analyze Lorentz violating thresholds in a 
rigorous manner. 



Rotation Invariance. Almost all work on thresholds to date has made the assumption that 
rotational invariance holds. If this invariance is broken, then our threshold theorems and results 
do not necessarily hold. For threshold discussions, we will assume that the underlying EFT is 
rotationally invariant and use the notation p = |7?|. 

Monotonicity. We will assume that the dispersion relations for all particles is monotonically 
increasing. This is the case for the mSME with small Lorentz violating coefficients if we work in 
a concordant frame. Mass dimension > 4 operators generate dispersion relations of the form 

E i = m i+ P * + f^P (88) 
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which do not satisfy this condition at momentum near the Planck scale if < 0. The turnover 
momentum pxo where the dispersion relation is no longer monotonically increasing is pro — 
(-2/(n/("))) 1 /("- 2 )£;p Z . The highest energy particles known to propagate ar the trans-GZK cos- 
mic rays with energy 10 -8 Epi- Hence unless 3> 1, pto is much higher than any relevant 
observational energy and we can make the assumption of monotonicity without loss of generality. 

High energy incoming particle. If there is a multi-particle in state, we will assume that 
one of the particles is much more energetic than all the others. This is the observational situation 
in reactions such as photon-photon scattering or pion production by cosmic rays scattering off the 
cosmic microwave background (the GZK reaction, see section l?i.5.6JI . 

6.5.3 Threshold theorems 

Eventually, any threshold analysis must solve for the threshold energy of a particular reaction. To 
do this, we must first know the appropriate kinematic configuration that applies at a threshold. 
Of use will be a set of threshold theorems that hold in the presence of Lorentz violation, which we 
state below. Variations on these theorems were derived in [BE] for single particle decays with n=2 
type dispersion and |215) for two in-two out particle interactions with general dispersion. Here we 
state the more general versions. 

Theorem 1: The configuration at a threshold for a particle with momentum p\ is the mini- 
mum energy configuration of all other particles that conserves momentum. 

Theorem 2: At a threshold all outgoing momenta are parallel to p\ and all other incoming 
momentum are anti-parallel. 



6.5.4 New threshold phenomena 
Asymmetric thresholds 

Asymmetric thresholds are thresholds where two outgoing particles with equal masses have un- 
equal momenta. This cannot occur in Lorentz invariant reactions. Asymmetric thresholds occur 
because the minimum energy configuration is not necessarily the symmetric configuration. To see 
this, let us analyze photon decay, where we have one incoming photon with momentum pi n and 
an electron/positron pair with momenta q\,q2- We will assume our Lorentz violating coefficients 
are such that the electron and positron have identical dispersion. 20 

Imagine that the dispersion coefficients for the electron and positron are negative and such 
that the electron/positron dispersion is given by the solid curve in Figure [5] We define the energy 
Esymm to be the energy when both particles have the same momentum q\ = qi = Pi n /2. This is not 
the minimum energy configuration, however, if the curvature of the dispersion relation (d 2 E/dp 2 ) 
at Pi n /2 is negative. If we add a momentum Aq to <?2 and — Ag to gi then we change the total 
energy by AE — AE2 — AE\. Since the curvature is negative, AEi > AE2 and therefore AE < 0. 
The symmetric configuration is not the minimum energy configuration and is not the appropriate 
configuration to use for a threshold analysis for all pi n . Note that part of the dispersion curve in 
Figure has positive curvature, as must be the case if at low energies we have the usual Lorentz 
invariant massive particle dispersion. If we were considering the constraints derivable when Pi n /2 
is small and in the positive curvature region, then the symmetric configuration would be the ap- 

20 This can easily be done by considering a CPT invariant EFT or choosing the appropriate helicity 
states/coefficients for a CPT violating EFT. 
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Figure 2: Total outgoing particle energy in symmetric and asymmetric configurations. 
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plicablc one. In general when it is appropriate to use asymmetric thresholds or symmetric ones 
depends heavily on the algebraic form of the outgoing particle Lorentz violation and the energy 
that the threshold must be above. The only general statement that can be made is that asymmetric 
thresholds are not relevant when the outgoing particles have n — 2 type dispersion modifications 
(either positive or negative) or for strictly positive coefficients at any n. For further examples of 
the intricacies of asymmetric thresholds, see [1541 1167] . 

Hard Cerenkov thresholds 



Related to the existence of asymmetric thresholds is the hard Cerenkov threshold, which also 
occurs only when n > 2 with negative coefficients. However, in this case both the outgoing and 
incoming particles must have negative coefficients. To illustrate the hard Cerenkov threshold, we 
consider photon emission from a high energy electron, which is the rotated diagram of the photon 
decay reaction. In Lorentz invariant physics, electrons emit soft Cerenkov radiation when their 
group velocity dE/dp exceeds the phase velocity cu/k of the electromagnetic vacuum modes in a 
medium. This type of Cerenkov emission also occurs in Lorentz violating physics when the group 
velocity of the electrons exceeds the low energy speed of light in vacuum. The velocity condition 
does not apply to hard Cerenkov emission, however, so to understand the difference we need to 
describe both types in terms of energy-momentum conservation. 

Let us quickly remind ourselves where the velocity condition comes from. The energy conser- 
vation equation (imposing momentum conservation) can be written as 

tj(k) = E(p)-E{p-k). (89) 

Dividing both sides by k and taking the soft photon limit k — > we have 

,. Lu(k) E(p)-E(p-k) dE 

lim — = hm — — r-^- = — . 90 

fe^o k fe^o k op 

Equation (|90|l makes clear that the velocity condition is only applicable for soft photon emission. 
Hard photon emission can occur even when the velocity condition is never satisfied if the photon 
energy-momentum vector is spacelike with n > 2 dispersion. As an example, consider an unmod- 
ified electron and a photon dispersion of the form to 2 — k 2 — k 3 /Epi. The energy conservation 
equation in the threshold configuration is 

m 2 m 2 k 2 .„„, 

P+ 2p- =p - k+ 2( p -k) +k '^EV l (91) 

where p is the incoming electron momentum. Introducing the variable x = k/p and rearranging 
we have 



vnfEpi 



= x(l-x). (92) 



Since all particles are parallel at threshold, x must be between and 1. The maximum value of the 
right hand side is 1/4, and so we see that we can solve the conservation equation if p > (Am^Epi) 1 ^ 3 , 
which is approximately 23 TeV. At threshold, x = 1/2 so this corresponds to emission of a hard 
photon with an energy of 11.5 TeV. 



Upper thresholds 



Upper thresholds do not occur in Lorentz invariant physics. It is easy to see that they are possible 
with Lorentz violation, however. In figure|31the region R in energy space spanned by Eo u t{Xk,Pi) 
is bounded below, since each individual dispersion relation is bounded below. However, if one can 
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adjust the dispersion E\(p{) freely, as would be the case if the incoming particle was a unique 
species in the reaction, then one can choose Lorentz violating coefficients such that E\{p\) moves 
in and out of R. 




Figure 3: An example of an upper and lower threshold. R is the region spanned by all Xj~ and 
E\{pi) is the energy of the incoming particle. Where E\{p{) enters and leaves R are lower and 
upper thresholds, respectively. 

As a concrete example consider photon decay, 7 — ► e++e~, with unmodified photon dispersion 
and an electron/positron dispersion relation of 



£ 2 =p 2 + m 2 -^p 2 + -^- (93) 
2E% 3 e pi 



chosen strictly for algebraic convenience. This dispersion relation has positive curvature every- 
where implying that the the electron and positron have equal momenta at threshold. The energy 
conservation equation, where the photon has momentum k is then 

.fc m 2 7to 2/3 , k 3 N 

t=2< 2 + — ^'"W (94) 
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which reduces to 

9 7m 2 / 3 9 k 3 , „ 

8to 2 — k 2 + — = 0. 95 

Ej* Epi 

<|95[) has two positive real roots, at k = (4 ± 2y/2)m 2 / 3 Ep{ 3 , corresponding to a lower and upper 
threshold at 14 and 82 TeV, respectively. Such a threshold structure would produce a deficit in the 
observed photon spectrum in this energy band. 21 Very little currently exists in the literature on 
the observational possibilities of upper thresholds. A complicated lower/upper threshold structure 
has been applied to the trans-GZK cosmic ray events, with the lower threshold mimicking the 
GZK-cutoff at 5 • 10 19 GeV and the upper entering below the highest events at 3 • 10 20 GeV [TBI] . 
The region of parameter space where such a scenario might happen is extremely small, however. 

Helicity decay 

In previous work on the Cerenkov effect based on EFT it has been assumed that left and right 
handed fermions have the same dispersion. As we have seen, however, this need not be the case. 
When the fermion dispersion is helicity dependent the phenomenon of helicity decay occurs. One 
of the helicities is unstable and will decay into the other as a particle propagates, emitting some 
sort of radiation depending on the exact process considered. Helicity decay has no threshold in 
the traditional sense; the reaction happens at all energies. However, below a certain energy the 
phase space is highly suppressed, so we have an effective threshold that practically speaking is 
indistinguishable from a real threshold. As an example, consider the reaction — > ep + 7, with 
an unmodified photon dispersion and the electron dispersion relation 

E 2 R = m 2 + p 2 (96) 
E 2 = m 2 +p 2 + f^p 2 (97) 

for right and left-handed electrons. Furthermore, assume that ff£ > 0. The opening up of the 
phase space can be seen by looking at the minimum and maximum values of the longitudinal 
photon momentum. The energy conservation equation is 

2 f(4) 2 

777 f 711 

p+ _ + ^ p . p - t + _ + | fi | (98 ) 

where p is the incoming momentum and k is the outgoing photon momentum. We have assumed 
that the transverse momentum is zero, which gives us the minimum and maximum values of k. k 
is assumed to be less than p\ one can check a posteriori that this assumption is valid. It can be 
negative, however, which is different from a threshold calculation where all momenta are necessarily 
parallel. Solving 1|98[) for fc m m, k ma x yields to lowest order in to, 

f(4) 

U - JeLP (gg) 



f (4) 2 

b - n JeLP 

m 2 + fe L 'p 2 

From (|99(l it is clear that when p 2 << m 2 /fj^ the phase space is highly suppressed, while for 
p 2 << to 2 / the phase space in k becomes of order p. The momentum p t h — (m 2 / f^) 1 ^ 2 acts 

21 Currently, the observed photon spectrum only extends to 50 TeV. Hence this type of lower/upper threshold 
structure is observationally indistinguishable from a simple lower threshold with current data. 
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as an effective threshoid, where the reaction is strongly suppressed below this energy. Constraints 
from helicity decay in the current literature |155) are complicated and not particularly useful. 
Hence we shall not describe them here, instead focussing our attention on the strict Cerenkov 
effect when the incoming and outgoing particle have the same helicity. For an in depth discussion 
of helicity decay constraints see [155) . 

6.5.5 Threshold constraints in QED 

With the general phenomenology of thresholds in hand, we now turn to the actual observational 
constraints from threshold reactions in Lorentz violating QED. We will continue to work in a 
rotationally invariant setting. Only the briefest listing of the constraints is provided here; for a 
more detailed analysis see (1541 11581 1155) . Most constraints in the literature have been placed by 
demanding that the threshold for an unwanted reaction is above some observed particle energy. 
As mentioned previously, a necessary step in this analysis is to show that the travel times of 
the observed particles are much longer than the reaction time above threshold. A calculation 
of this for the vacuum Cerenkov has been done for QED with dimension four Lorentz violating 
operators in 224 . More generally, a simple calculation shows that the energy loss rate above 
threshold from the vacuum Cerenkov effect rapidly begins to scale as e 2 AE n / Ep\ 2 where A is a 
coefficient that depends on the coefficients of the Lorentz violating terms in the EFT. Similarly, 
the photon decay rate is e 2 AE n ~ 1 / EpJ 2 . In both cases the reaction times for high energy particles 
are roughly (e 2 A) -1 EpJ 2 / E n ~ x , which is far shorter than the required lifetimes for electrons and 
photons in astrophysical systems for n = 2, 3. 22 The lifetime of a high energy particle in QED 
above threshold is therefore short enough that we can establish constraints simply by looking at 
threshold conditions. 

Photon decay Lorentz violating terms can be chosen such that photons become unstable to 
decay into electron-positron pairs |152| . We observe 50 TeV photons from the Crab nebula. There 
must exist then at least one stable photon polarization. The thresholds for n = 2, 3 dispersion have 
been calculated in [154] . Demanding that these thresholds are above 50 TeV yields the following 
best constraints. 

For n = 2 with CPT preserved we have ffl - / e (2) < 4m 2 /p 2 h = 4 • 10" 16 If we set d = 

in ()40|) so that there is no helicity dependence, this translates to the constraint fcj?/2 + c < 4- 10~ 16 . 
If d then both helicities of electrons/positrons must satisfy this bound since the photon has 
a decay channel into every possible combination of electron/positron helicity. The corresponding 
limit is k F /2+ {c±d) < 4- 10~ 16 . 

For n = 3 the situation is a little more complicated, as we must deal with photon and electron 
helicity dependence, positron dispersion, and the possibility of asymmetric thresholds. The 50 
TeV Crab photon polarizations are unknown, so only the region of parameter space in which both 
polarizations decay can be excluded. We can simplify the problem dramatically by noting that the 
birefringence constraint on £ in 1)42)1 is |£| < 10~ 4 |127) . The level of constraints from threshold 
reactions at 50 TeV is around 10 -2 [15211167) . Since the birefringence constraint is so much stronger 
than threshold constraints, we can effectively set £ = and derive the photon decay constraint in 
the region allowed by birefringence. 

With this assumption, we can derive a strong constraint on both rjn and tjl by considering 
the individual decay channels 7 — > eJp + and 7 — > + ep, where L,R stand for the helicity. 
For brevity, we shall concentrate on 7 — > e~p + e£, the other choice is similar. The choice of a 
right-handed electron and left-handed positron implies that both particle's dispersion relations are 

22 For n = 4 no QED particles reach energies high enough to provide constraints. The only particles of the required 
energy are ultra-high energy cosmic rays or neutrinos. Assuming the cosmic rays are protons, the corresponding 
reaction time for Cerenkov emission is 10 — 17 sec. 
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functions of only /' ; „ and hence rju (see Sec. I4.1.4|l . The matrix element can be shown to be large 



leR 

enough for this combination of helicities that constraints can be derived by simply looking at the 
threshold. Imposing the threshold configuration and momentum conservation and substituting in 
the appropriate dispersion relations the energy conservation equation becomes 

m 2 f (3) m 2 f (3) 

h =* + % + kr/ +k -* + W^j-k-, ik -<' ) (100) 

where k is the incoming photon momentum and p is the outgoing electron momentum. Cancelling 
the lowest order terms and introducing the variable z = 2p/k — 1 this can be rewritten as 

- 4m2Epi (101) 



To find the minimum energy configuration we must minimize the right hand side of i|l(Jlfl with 
respect to z (keeping the right hand side positive). We note that since the range of z is between 
-1 and 1, the right hand side of (|101|) can be positive for both positive and negative fjj}, which 
implies that the bound will be two sided. 

As an aside, it may seem odd that photon decay happens at all when the outgoing particles 
have opposite dispersion modifications since the net effect on the total outgoing energy might 
seem to cancel. However, this is only the case if both particles have the same momenta. We can 
always choose to place more of the incoming momentum into the outgoing particle with a negative 
coefficient, thereby allowing the process to occur. This reasoning also explains why the bound 
is two sided, as the threshold configuration gives more momentum to whichever particle has a 
negative coefficient. 

Returning to the calculation of the threshold, minimizing (|101|) . we find that the threshold 
momentum is 

m 2 Epi , 



^thresh 



A3) 
J e 



(102) 
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The absolute value here appears because we find the minimum positive value of (|101|) . Placing 

kthresh at 50 TeV yields the constraint \f% !\ < 0.25 and hence \t)r\ < 0.125. The same procedure 
applies in the opposite choice of outgoing particle helicity, so r\p obeys this bound as well. 

Vacuum Cerenkov The 50 TeV photons observed from the Crab nebula are believed to be 
produced via inverse Compton (IC) scattering of charged particles off the ambient soft photon 
background. 23 If one further assumes that the charged particles are electrons, it can then be 
inferred that 50 TeV electrons must propagate. However, only one of the electron helicities may 
be propagating, so we can only constrain one of the helicities. 

For n = 2 the constraint is — f^p < m 2 /p 2 h — 10~ 16 |154) . where /i 2 "* is the coefficient for 
one of the electron helicities. In terms of the mSME parameters this condition can be written as 
c — \d\ — kp/2 < 10~ 16 . For n = 3 an added complication arises. If we consider just electrons as 

(3) (3) 

the source of the 50 TeV photons, then we have that either or /\j must satisfy 
a) /i 3 > < 0.012 for / e (3) > and ^ 3) > -3/i 3) , (103) 



b) 4 3) > / p (3) - 0.048 - 2V(0.024) 2 - 0.048/ e (3) for f\ 3) < -3/ e (3) < or /^ 3) < / e (3) < 



23 The theoretical model for the non-thermal emission of the Crab nebula is called the synchrotron-self-Compton 
(SSC) model. For recent fits of data to the SSC spectrum see 
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and the translation to £ and t]r.l is as before. Note that for the range of £ allowed by birefringence, 
the relevant constraint is r\B. < 0.012 or r)L < 0.012. 

A major difficulty with the above constraint is that positrons may also be producing some of 
the 50 TeV photons from the Crab nebula. Since positrons have opposite dispersion coefficients in 
the n = 3 case, there is always a charged particle able to satisfy the Cerenkov constraint. Hence 
by itself, this IC Cerenkov constraint can always be satisfied in the Crab and gives no limits at all. 
However, as we shall see in section the vacuum Cerenkov constraint can be combined with the 
synchrotron constraint to give an actual two-sided bound. 

Photon annihilation The high energy photon spectrum (above 10 TeV) from astrophysical 
sources such as Markarian 501 and 421 has been observed to show signs of absorption due to 
scattering off the IR background. While this process occurs in Lorentz invariant physics, the 
amount of absorption is affected by Lorentz violation. The resulting constraint is not nearly as 
clear cut as in the photon decay and Cerenkov cases, as the spectrum of the background IR photons 
and the source spectrum are both important, neither of which is entirely known. Various authors 
have argued for different constraints on the n — 3 dispersion relation, based upon how far the 
threshold can move in the IR background. The constraints vary from 0(1) to O(10). However, 
none of the analyses take into account the EFT requirement for n = 3 that opposite photon 
polarization have opposite Lorentz violating terms. Such an effect would cause one polarization 
to be absorbed more strongly than in the Lorentz invariant case and the other polarization to be 
absorbed less strongly. The net result of such a situation is currently unknown, although current 
data from blazars suggest that both polarizations must be absorbed to some degree |263) . Since 
even at best the constraint is not competitive with other constraints, and since there is so much 
uncertainty about the situation, we will not treat this constraint in any more detail. For discussions 
see [TCTII7|. 

6.5.6 The GZK Cutoff and Ultra-high energy cosmic rays 

The GZK Cutoff Ultra-high energy cosmic rays (UHECR), if they are protons, will interact 
strongly with the cosmic microwave background and produce pions, p + j — * p + ir°, losing energy 
in the process. As the energy of a proton increases, the GZK reaction can happen with lower and 
lower energy CMBR photons. At very high energies (5- 10 19 eV), the interaction length (a function 
of the power spectrum of interacting background photons coupled with the reaction cross section) 
becomes of order 50 Mpc. Since cosmic ray sources are probably at further distances than this, 
the spectrum of high energy protons should show a cutoff around 5 • 10 19 eV 135 281 . A number 
of experiments have looked for the GZK cutoff, with conflicting results. AGASA found trans-GZK 
events inconsistent with the GZK cutoff at 2.5cr 96 while Hi-Res has found evidence for the GZK 
cutoff (although at a lower confidence level, for a discussion see 263]). New experiments such 
as AUGER [113] may resolve this issue in the next few years. Since Lorentz violation shifts the 
location of the GZK cutoff, significant information about Lorentz violation (even for n — 4 type 
dispersion) can be gleaned from the UHECR spectrum. If the cutoff is seen then Lorentz violation 
will be severely constrained, while if no cutoff or a shifted cutoff is seen then this might be a 
positive signal. 

For the purposes of this review, we will assume that the GZK cutoff has been observed and 
describe the constraints that follow. We can estimate their size by noting that in the Lorentz 
invariant case the conservation equation can be written as 

(p + kf = (m, + m v f (104) 

as the outgoing particles are at rest at threshold. Here p is the UHECR proton 4-momentum and 
k is the soft photon 4=momentum. At threshold the incoming particles are anti-parallel, which 
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gives a threshold energy for the GZK reaction of 



E GZK ^ 3 ■ l(r u GeV • 




2.7K 



) 



(105) 



where ojq is the energy of the CMBR photon. The actual GZK cutoff occurs at 5- 10 19 eV due to the 
tail of the CMBR spectrum and the particular shape of the cross section (the A resonance). From 
this heuristic threshold analysis, however, it is clear that Lorentz violation can become important 
when the modification to the dispersion relation is of the same order of magnitude as the proton 
mass. For n = 2 dispersion, a constraint of /i- — < O(10~ 23 ) was derived in |88l 1871 ITU] , 
The case of n = 3 dispersion with f {3) = ff ] was studied in [TSUI HS2 EHU 1351 1371 1391 1251 D51 
1231 Cnni E21 CSZl El, while the possibility of fi 3A) ^ was studied in [T53]. A simple 

constraint |154) can be summarized as follows. If we demand that the GZK cutoff is between 
2 • 10 19 eV and 7 • 10 19 eV then for f {3) = ff ] we have |/^ 3) | < O(10" 14 ). If fi 3) ^ /j 3) then there 
is a wedge shaped region in the parameter space that is allowed |154) . 

The numerical values of these constraints should not be taken too literally. While the order of 
magnitude is correct, simply moving the value of the threshold for the proton that interacts with a 
CMBR photon at some energy does not give accurate numbers. GZK protons can interact with any 
photon in the CMBR distribution above a certain energy. Modifying the threshold modifies the 
phase space for a reaction with all these photons in the region to varying degrees, which must be 
folded in to the overall reaction rate. Before truly accurate constraints can be calculated from the 
GZK cutoff a more detailed analysis to recompute the rate in a Lorentz violating EFT considering 
the particulars of the background photon distribution and A-resonance must be done. However, 
the order of magnitude of the constraints above is roughly correct. Since they are so strong, the 
actual numeric coefficient is not particularly important. 2 

Another difficulty with constraints using the GZK cutoff is the assumption that the source 
spectrum follows the same power law distribution as at lower energies. It may seem that proposing a 
deviation from the power law source spectrum at that energy would be a conspiracy and considered 
unlikely However, this is not quite correct. A constraint on /' n J will, by the arguments above, be 
such that the Lorentz violating terms are important only near the GZK energy - below this energy 
we have the usual Lorentz invariant physics. However, such new terms could then also strongly 
affect the source spectrum only near the GZK energy. Hence the GZK cutoff could vanish or be 
shifted due to source effects as well. Unfortunately, we have little idea as to the mechanism that 
generates the highest energy cosmic rays, so we cannot say how Lorentz violation might affect their 
generation. In summary, while constraints from the position of the GZK cutoff are impressive and 
useful, their actual values should be taken with a grain of salt since a number of unaccounted for 
effects may be tangled up in the GZK cutoff. 

UHECR Cerenkov A complimentary constraint to the GZK analysis can be derived by recog- 
nizing that 10 19 — 10 20 eV protons reach us - a vacuum Cerenkov effect must be forbidden up to the 
highest observed UHECR energy |88l 11541 1119) . The direct limits from photon emission, treating 



a 5 • 10 19 eV proton as a single constituent are ff ] - f!f ] < 4 • 10" 22 HHU EH1 ESI for n = 2, 25 
fv 3) ~ A 3) < 0(1CT 14 ) for n = 3 [T53], and - f\ 4} < 0(1CT 5 ) for n = 4 [T53]. Equivalent 
bounds on Lorentz violation in a conjectured low energy limit of loop quantum gravity have also 
been derived using UHECR Cerenkov |190) . 



This would change, of course, if a shift in the cutoff was recorded, as then the detailed dependence of the cutoff 
location on the Lorentz violating physics would be important. 

25 The bounds quoted in I5S1ITT^1 vary from this bound by roughly an order of magnitude, as the proton energy 
is taken to be 10 20 eV rather than right at the GZK cutoff. 
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Cerenkov emission for UHECR has been used most extensively in |119j . where two-sided limits 
on Lorentz violating dimension 4,5, and 6 operators for a number of particles are derived. The 
argument is as follows. If we view a UHECR proton as actually a collection of constituent partons 
(i.e. quarks, gauge fields, etc.) then the dispersion correction should be a function of the corrections 
for the component partons. By evaluating the parton distribution function 26 for protons and other 
particles at high energies, one can get two sided bounds by considering multiple reactions, in the 
same way one obtains two sided bounds in QED. As a simple example, consider only dimension 
four rotationally invariant operators (i.e. n = 2 dispersion) and assume that all bosons propagate 
with speed one while all fermions have a maximum speed of 1 — e. Let us take the case e < 0. A 
proton is about half fermion and half gauge boson, while a photon is 80 percent gauge boson and 20 
percent fermion. The net effect, therefore, is that a proton travels faster than a photon and hence 
Cerenkov radiates. Demanding that a 10 20 eV proton not radiate yields the bound e > — 10~ 23 , 
similar to the standard Cerenkov bound above. 

If instead e > 0, then e + e~ pair emission becomes possible as electrons and positrons are 85 
percent fermion and 15 percent gauge boson. Pair emission would also reduce the UHECR energy, 
so one can demand that this reaction is forbidden as well. This yields the bound e < 10~ 23 . 
Combined with the above bound we have |e| < 10~ 23 , which is a strong two sided bound. The 
parton approach yields two-sided bounds on dimension six operators of order I/*- 4 -*! < O(10~ 2 ) 
for all constituent particles, depending on the assumptions made about equal parton dispersion 
corrections. Bounds on the coefficients of CPT violating dimension five operators are of the order 
10 -15 . 27 For the exact constraints and assumptions, see |119) . Note that if one treated electrons, 
positrons, and protons as the fundamental constituents with only n = 2 dispersion and assigned 
each a common speed 1 — e one would obtain no constraints. Therefore the parton model is more 
powerful. However, for higher dimension operators that yield energy dependent dispersion, simply 
assigning electrons and protons equal coefficients f( n > does yield comparable constraints. Finally, 
we comment that |119) does not explicitly include possible effects such as SUSY that would change 
the parton distribution functions at high energy. 

6.5.7 Gravitational Cerenkov 

High energy particles travelling faster than the speed of graviton modes will emit graviton Cerenkov 
radiation. The authors of |224) have analyzed the emission of gravitons from a high energy particle 
with n — 2 type dispersion and find the rate to be 

dE Gp 4 , n9 , , 

^^(c P -l) 2 (106) 

where c p is the speed of the particle and G is Newton's constant. We have normalized the speed of 
gravity to be one. The corresponding constraint from the observation of high energy cosmic rays 
is c p — 1 < 2 ■ 10~ 15 . This bound assumes that the cosmic rays are protons, uses the highest record 
energy 3 • 10 20 eV, and assumes that the protons have travelled over at least 10 kpc. Furthermore, 
the bound assumes that all the cosmic ray protons travel at the same velocity, which is not the 
case if CPT is violated or d ^ in the mSME. 

The corresponding bounds for n = 3,4 type dispersion are not known, but one can easily 
estimate their size. The particle speed is approximately 1 + f^ n ' (E j 'Epi) n ~ 2 . For a proton at 
an energy of 10 20 eV (10~ 8 Epi) the constraint on the coefficient is then of O(10 -7 ). Note 

26 This approach requires that one can evaluate the parton distribution functions up to UHECR energies with- 
out errors from the uncertainty in low energy parton distribution functions or effects from new physics becoming 
appreciable. We caution that such effects can adjust the constraints, as the authors of 11191 note. 

27 Dimcnsion five CPT violating operators yield helicity dependent dispersion. To derive bounds on these opera- 
tors rrtTn assumes that the coefficients are roughly equal. 
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though, that in this case only one of the UHECR protons must satisfy this bound due to helicity 
dependence. Similarly, the n = 4 bound is of O(10). 

Equation (|106|) only considers the effects of Lorentz violation in the matter sector which give 
rise to a difference in speeds, neglecting the effect of Lorentz violation in the gravitational sector. 
Specifically, the analysis couples matter only to the two standard graviton polarizations. However, 
as we shall see in Sec. 17.11 consistent Lorentz violation with gravity can introduce new gravitational 
polarizations with different speeds. In the aether theory ( section 14. 4JI there are three new modes, 
corresponding to the three new degrees of freedom introduced by the constrained aether vector. 
The corresponding Cerenkov constraint from possible emission of these new modes has recently 
been analyzed in I105j . Demanding that high energy cosmic rays not emit these extra modes and 
assuming no significant Lorentz violation for cosmic rays yields the bounds 

- ci - c 3 < i x lcr 15 



(c 1 +c 3 ) 2 (c 2 + 3c 1 c 3 -2c 4 ) 1 4 x 1Q _ ai 



(c 3 -c 4 ) 2 ^ 1 w 1A _ 30 



I 
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< 1 x 10 



\d + C 4 | 

C4 ' C2 ' C3 < 3 x 10- 19 (107) 

Cl 

on the coefficents in (I50|) . The next to last bound requires that (c 4 — C2 — C3)/ci > 10~ 22 . If, as the 
authors of [105] argue, no gravity-aether mode can be superluminal, then these bounds imply that 
every coefficient is generically bounded by |c$| < 10 -15 . There is, however, a special case given by 
C3 = — ci, C4 = 0, C2 = ci/(l — 2ci) where all the modes propagate at exactly the speed of light and 
hence avoid this bound. 



6.6 Threshold reactions in other models 
6.6.1 Thresholds and DSR 

Doubly special relativity modifies not only particle dispersion relation but also the form of the 
energy conservation equations. The situation is therefore very different from that in EFT. The 
first difference between DSR and EFT is that DSR evades all of the photon decay and vacuum 
Cerenkov constraints that give strong limits on EFT Lorentz violation. Since there is no EFT type 
description of particles and fields in a DSR framework, one has no dynamics and cannot calculate 
reaction rates. However, one still can use the DSR conservation laws to analyze the threshold 
kinematics. By using the pseudo-momentum and and energy tt, e one can show that if a reaction 
does not occur in ordinary Lorentz invariant physics, it does not occur in DSR |146j . Physically, 
this is obvious. If the vacuum Cerenkov effect for, say. electrons began to occur at some energy E th , 
in a different reference frame the reaction would occur at some other energy EL, as the threshold 
energy is not an invariant. Therefore frames could be distinguished by labelling them according 
to the energy when the vacuum Cerenkov effect for electrons begins to occur. This violates the 
equivalence of all inertial frames that is postulated in DSR theories. 

A signal of DSR in threshold reactions would be a shift of the threshold energies for reactions 
that do occur, such as the GZK reaction or 7-ray annihilation off the infrared background |21) . 
However, the actual shift of threshold energies due to DSR is negligible at the level of sensitivity 
we have with astrophysical observations 21 . Hence DSR cannot be ruled out or confirmed by any 
threshold type analysis we currently have. The observational signature of DSR would therefore be 
a possible energy dependence of the speed of light fsection l6~2*)l without any appreciable change in 
particle thresholds |258| . 
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6.6.2 Thresholds and non-systematic dispersion 

Similar to DSR, the lack of dynamics in the non-systematic dispersion framework of Sec. 13.51 
makes it more problematic to set bounds on the parameters f( n \ In [1601 1121 1111 123) . the authors 
assume that the net effect of spacetime foam can be derived by considering energy conservation 
and non-systematic dispersions at a point. There is a difficulty with this, which we shall address, 
but for now let us assume that this approach is correct. 

As an example of the consequences of non-systematic dispersion let us consider the analysis 
of the GZK reaction in jllj . The authors consider n = 3 nonsystematic dispersion relations with 
normally distributed coefficients fp t 4 that can take either sign and have a variance of 0(1). Looking 
solely at the kinematical threshold condition, they find that all cosmic ray protons would undergo 
photo-pion production at energies above 10 15 eV. This is perhaps expected, as the energy scale 
that an n = 3 term becomes important is E « (m 2 E P i/ /( 3 ))V3 ~ 10 15 e y for /( 3 ) of 0(1). There 

(3) 

is a large region of the fp^ parameter space that is susceptible to the vacuum Cerenkov effect with 
pion emission |154| and hence a significant amount of the time the random coefficients will fall in 
this region of parameter space. If an ultra-high energy proton can emit a pion without scattering 
off of the CMBR, then certainly it can scatter as well, which implies that GZK reaction is also 
accessible. This same type of argument can be rapidly extended to n — 4 dispersion, yielding a 
cutoff in the spectrum at 10 18 eV. The n = 4 cutoff could easily be pushed above GZK energies if 
the coefficients had a variance slightly less than 0(1). In short, since we see high energy cosmic 
rays at energies of 10 20 eV, the results of [HH E2 El E3 imply that we could not have n = 3 
non-systematic dispersion unless the coefficients are <C O(l), while for n — 4 the coefficients would 
only have to be an order of magnitude or two below 0(1). 

We now return to a possible problem with this type of analysis, which has been raised in |43| . 
Performing threshold analyses on non-systematic dispersion assumes that energy-momentum con- 
servation can be applied with a single fluctuation (i.e. the reaction effectively happens at a point). 
It further assumes that the matrix element is roughly unchanged. In GZK or Cerenkov reactions, 
however, one of the outgoing particles is much softer than the incoming particle. In this situation 
the interaction region is much larger than the de Broglie wavelength of the high energy incoming 
particle, which means that many dispersion fluctuations will occur during the interaction. The 
amplitude of low energy emission in regular quantum field theory changes dramatically in this sit- 
uation (e.g. Bremsstrahlung with a rapidly wiggling source) as opposed to the case in which there 
is only one fluctuation (e.g. the Cerenkov effect). The above approach, modified conservation plus 
unchanged matrix element/rate when the reaction is allowed, is not correct when a low energy 
particle is involved. If the outgoing particle has an energy comparable to the incoming particle, 
then it may be possible to avoid this problem. However, in this case the reverse reaction is also 
kinematically possible with a different fluctuation of the same order of magnitude, so it is unclear 
what the net effect on the spectrum should be. Note, finally, that these arguments only concern 
the rate of decay - the conclusion that high energy particles would decay in this framework is 
unchanged. 

6.7 Synchrotron radiation 

Ref. |153) noticed that the synchrotron emission from the Crab nebula (and other astrophysical 
objects) is very sensitive to n = 3 modified dispersion relations. The logic behind the constraint 
in |153| is as follows. The Crab nebula emits electromagnetic radiation from radio to multi-TeV 
frequencies. As noted in section l6.5.5l the spectrum is well fit over almost the entire frequency range 
by the synchrotron self-Compton model. If this model is correct, as seems highly probable, then 
the observed radiation from the Crab at 100 MeV is due to synchrotron emission from high energy 
electrons and/or positrons. The authors of [153] . in the context of effective field theory, argued 
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that the maximum frequency to c of synchrotron radiation in the presence of Lorentz violation is 
given by 

3 eB-f 3 

L ° c = ~2~E~ (108) 
where e is the charge, B is the magnetic field, and 7, E are the gamma factor and energy of the 
source particle. The derivation of l|lU8fl was challenged in [7S]. More detailed calculations 222 
111111221] show that in the case of the Crab nebula (|108|) is correct, although the argument of |153) 
does not necessarily hold in general 28 . Assuming the source particles are electrons, if < then 
there is a maximum electron velocity for any energy particle and hence a maximum possible value 
for lo c . The maximum frequency must be above 100 MeV in the Crab, which leads to a constraint 
of either f^' R or 

fe 3 l > -7 • 10~ 8 , i.e. at least one of the electron parameters must be above this 

value. 

The analysis of |153| does not take into account the possibility that the high energy synchrotron 
emission could be due to positrons, which may also be generated near the pulsar. This is an im- 
portant possibility, since in the EFT that gives rise to J*- 3 -* terms for electrons, (I42|) . the positron 
has an opposite dispersion modification. Hence there is always some charged particle in the Crab 
with a dispersion modification that evades the synchrotron constraint. The possibility that there 
are two different populations, one electrons and one positrons, that contribute to the overall spec- 
trum would be a departure from the synchrotron-self Compton model, which presupposes only one 
population of particles injected into the nebula. However, such a possibility cannot be ruled out 
without more detailed modelling of the Crab nebula and a better understanding of how the initial 
injection spectrum of particles from the pulsar is produced. 29 

The possible importance of positrons in the Crab implies that the synchrotron constraint is 
always satisfied when considered by itself. However, the synchrotron constraint can be combined 
with Cerenkov constraints to create a two sided bound |155j . Essentially, the SSC model is such 
that whatever species of particle is producing the synchrotron spectrum must also be responsible 
for the inverse Compton spectrum. Hence at least one helicity of electron or positron must satisfy 
both the vacuum Cerenkov and synchrotron constraints. This is not automatically satisfied in 
EFT and constitutes a true constraint. For more discussion see |155) . The combined synchrotron, 
threshold, birefringence, and time of flight constraints are displayed in Figure^] (taken from |155| ^l . 
In this plot r)± is equal to f^ R L ^ (not to be confused with the "qn^L in H42JI which differs by a 
factor of two). 

Finally, note that if the effective field theory is CPT conserving then positrons/electrons have 
the same dispersion relation. So, for n — 2 dispersion the 100 MeV synchrotron radiation from 
the Crab yields a parallel constraint of roughly / + > — 10~ 22 for at least one helicity of elec- 
tron/positron. 



6.8 Neutrinos 

Neutrinos can provide excellent probes of Lorentz violation, as their mass is much smaller than any 
other known particle. To see this consider the modified dispersion framework. For an electron with 
n = 3, 4 dispersion the energies at which Lorentz violation can become appreciable are at 10, 10 5 
TeV respectively. However, for a neutrino with mass even at 1 eV the corresponding energies 
are only 1 GeV for n — 3 and 1 TeV for n — 4, well within the realm of accelerator physics. 
The most sensitive tests of Lorentz violation in the neutrino sector come from neutrino oscillation 

28 For a discussion of synchrotron radiation in other models with Planck suppressed dispersion corrections see 11281 . 

29 It has recently been suggested that Lorentz violation might play a role in the formation of the pulsar itself, 
specifically in large pulsar anomalous velocities 11921 . 
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Figure 4: Constraints on LV in QED at n = 3 on a log-log plot. For negative parameters minus the 
logarithm of the absolute value is plotted, and region of width 10~ 10 is excised around each axis. 
The constraints in solid lines apply to £ and both rj± , and are symmetric about both the £ and the 
X] axis. At least one of the two pairs (rf± , £) must lie within the union of the dashed bell-shaped 
region and its reflection about the £ axis. Intersecting lines are truncated where they cross. 
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experiments, which we now describe. For a more comprehensive overview of neutrino mixing, see 
for example |1U2II165] . 



6.8.1 Neutrino oscillations 

Lorentz violating effects in the neutrino sector have been considered by many authors [EE1 11821 
HHU HHOI E31 ESI HHI HnH E3 To illustrate how Lorentz violation affects neutrino propagation, 
we consider the simplest case, where the limiting speeds for mass eigenstates of the neutrino are 
different, i.e. neutrinos have dispersions 

E 2 = (l + f^)p 2 +m 2 (109) 

where i denotes the energy eigenstate. In this case, the energy eigenstates are also the mass 
eigenstates (this is not necessarily the case with general Lorentz violation). This is a special case 
of the neutrino sector of the mSME if we assume that c°° is flavor diagonal and is the only non-zero 
term. For relativistic neutrinos, we can expand the energy to be 

f (2) m 2 

Now consider a neutrino produced via a particle reaction in a definite flavor eigenstate / with 
energy E. We denote the amplitude for this neutrino to be in a particular energy eigenstate i by 
the matrix Un where J^i UjJJu = <5rj- The amplitude for the neutrino to be observed in another 
flavor eigenstate J at some distance L, T from the source is then 

= E U]^ l{ET - pL) Uu « £ U] % e-^ E+ ^ L U H (111) 

i i 

for relativistic neutrinos. If we define an "effective mass" Ni as 

= m\ + fifE 2 (112) 
then the probability Pjj = \Au\ 2 can be written as, 

Pu = Su - £ 4F LHjSl n 2 ( J + 2G IJijS in( -^g- J (113) 

i,j>i ^ ' ^ ' 

where SNij — Nf — Nj and F]jij, Guij are functions of the U matrices. 

We can immediately see from Ijll3ll that Lorentz violation can have a number of consequences 
for standard neutrino oscillation experiments. The first is simply that neutrino oscillation still 
occurs even if the mass is zero. In fact, some authors have proposed that Lorentz violation could 
be partly responsible for the observed oscillations [181| . Oscillations due to the type of Lorentz vi- 
olation above vary as EL |181| . Current data support neutrino oscillations that vary as a function 
of L/E jSSl 1125) so it seems unlikely that Lorentz violation could be the sole source of neutrino os- 
cillations. It is possible, however, that Lorentz violation may explain some of the current problems 
in neutrino physics by giving a contribution in addition to the mass term. For example it has been 
proposed in |182j that Lorentz violation might explain the LSND (Liquid Scintillator Neutrino 
Detector) anomaly |36|. 30 . which is an excess of — > events that cannot be reconciled with 
other neutrino experiments |129) . We note that the above model for Lorentz violating effects in 

30 There are other attempts to explain the LSND results with CPT violation 1231 13T1 . However, these CPT 
violating models are not directly correlated with Lorentz violation as they may involve non-local field theories. It 
is also seems likely that the LSND result is simply incorrect. 
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neutrino oscillations is perhaps the simplest case. In the neutrino sector of the mSME there can 
be more complicated energy dependence, directional dependence, and new oscillations that do not 
occur in the standard model. For a discussion of these various possibilities see |18(J) . 

The difference in speeds between electron and muon neutrinos was bounded in |88j to be 
l/y? — fvj}\ < 10~ 22 . Oscillation data from Super Kamiokande have improved this bound to 
O(10~ 24 ) |114) . Current neutrino oscillation experiments are projected to improve on this by three 
orders of magnitude, giving limits on maximal speed differences of order 10 -25 |12fi| . For compar- 
ison, the time of flight measurements from supernova 1987A constrain \ — fly \ < 10 -8 |265| . 
Neutrino oscillations are sensitive enough to directly probe non-renormalizable Lorentz violating 
terms. In |69| current neutrino oscillation experiments are shown to yield bounds on dimension five 
operators stringent enough that the energy scale suppressing the operator must be a few orders 
of magnitude above the Planck energy. Such operators are therefore very unlikely in the neutrino 
sector. Ultra-high energy neutrinos, when observed, will provide further information about neu- 
trino Lorentz violation. For example, flavor oscillations of ultra-high energy neutrinos at 10 21 eV 
propagating over cosmic distances would be able to probe Lorentz violating dispersion suppressed 
by seven powers of Epi anc k |83| (or more if the energies are even higher). 

Additionally, neutrino Lorentz violation can modify the energy thresholds for reactions involv- 
ing neutrinos, which can have consequences for the expected flux of ultra-high energy neutrinos for 
detectors such as ICECUBE. The expected flux of ultra-high energy neutrinos is bounded above by 
the Bahcall-Waxman bound |273j if the neutrinos are produced in active galactic nuclei or gamma 
ray bursters. It has been shown ^S] that Lorentz violation can in fact raise (or lower) this bound 
significantly. A higher than expected ultra-high energy neutrino flux therefore could be a signal of 
Lorentz violation. 

6.8.2 Neutrino Cerenkov effect 

Finally, neutrinos can also undergo a vacuum Cerenkov effect. Even though a neutrino is neutral 
there is a non-zero matrix element for interaction with a photon as well as a graviton. Graviton 
emission is very strongly suppressed and unlikely to give any useful constraints. The matrix 
element for photon emission, while small, is still larger than that for graviton emission and hence 
the photon Cerenkov effect is more promising. The photon-neutrino matrix element can be split 
into two channels, a charge radius term and a magnetic moment term. The charge radius interaction 
is suppressed by the W mass, leading to a reaction rate too low for current neutrino observatories 
such as AMANDA to constrain n — 3, 4 Lorentz violation. However, the rate from the charge 
radius interaction scales strongly with energy, and it has been estimated |154) that atmospheric 
PeV neutrinos may provide good constraints on n — 3 Lorentz violation. The magnetic moment 
interaction has not yet been conclusively analyzed, so possible constraints from the magnetic 
moment interaction are unknown. In Lorentz invariant physics, the magnetic moment term is 
suppressed by the small neutrino mass, so energy loss rates are likely small. However, it should 
be noted that some Lorentz violating terms in an effective field theory give rise to effective masses 
that scale with energy. These might be much larger than the usual neutrino mass at high energies, 
yielding a large neutrino magnetic moment. 

6.9 Phase coherence of light 

An interesting and less well known method of constraining non-systematic Lorentz violation is 
looking at Airy rings (interference fringes) from distant astrophysical objects. In order for an 
interference pattern from an astrophysical source to be observed, the photons reaching the detector 
must be in phase across the detector surface. However, if the dispersion relation is fluctuating, 
then the phase velocity = co/k is also changing. If the fluctuations are uncorrelated then an 
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initially in phase collections of photons will lose phase coherence as they propagate. Uncorrelated 
fluctuations are reasonable, since for most of their propagation time, photons that strike across a 
telescope mirror are separated by macroscopic distances. Observation of Airy rings implies that 
photons are in phase and hence limits the fluctuations in the dispersion relation [2041 12451 1233 . 
The aggregate phase fluctuation is given by p33 

t n—1 ]~)3—n 

Acp = 2nf^^-^- (114) 

A 

where Lpi is the Planck length, D is the distance to the source, and A is the wavelength of the 
observed light. This technique was originally applied by |2U4) . but the magnitude of the aggregate 
phase shift was overestimated. PKS1413+135, a galaxy at a distance of 1.2 Gpc, shows Airy rings 
at a wavelength of 1.6/im. Demanding that the overall phase shift is less than 2tt yields 0(1) 
constraints for n — 5/2 and constraints of order 10 9 for n = 8/3. Hence this type of constraint 
is only able to minimally constrain Lorentz violating non-systematic models. In principle the 
frequency of light used for the measurement can be increased, however, in which case this type 
of constraint will improve. However, Ref. has argued that other effects mask the loss of phase 
coherence from quantum gravity, making even this approach uncertain. 
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7 Gravitational observations 



So far we have restricted ourselves to Lorentz violating tests involving matter fields. It is also 
possible that Lorentz violation might manifest itself in the gravitational sector. There are three 
obvious areas where the consequences of such Lorentz violation might manifest itself: gravitational 
waves, cosmology, and post-Newtonian corrections to weak field general relativity. 

7.1 Gravitational waves 

In the presence of dynamical Lorentz violation, where the entire action is diffcomorphism invariant, 
one generically expects new gravitational wave polarizations 31 . The reason is simple. Any dynami- 
cal Lorentz violating tensor field must have kinetic terms involving derivatives of the form V /i J7 Q ' 3 " 
where U a @ "' is the Lorentz violating tensor. Furthermore, U must take a non-zero vacuum expec- 
tation value if it violates Lorentz invariance. At linear order in the perturbations h a p, u a P--- (where 
g a p = Va/3 + h a p, U a/3 --- =< U a/3 --- > +u Q ^ - ), the connection terms in the covariant derivative are 
also first order, for example d a hj3 1 < U^ S '" >. Upon varying the linearized metric, these terms 
contribute to the graviton equations of motion. The extra terms in the graviton equations give rise 
to new solutions. Since the potential that forces U to take a non-zero vacuum expectation value 
must involve the metric, variations in U are usually coupled to metric variations, implying that 
the new graviton modes mix with excitations of the Lorentz violating tensor fields. 

There is a large literature on gravitational wave polarizations in theories of gravity other than 
general relativity. For a thorough discussion, see |277j and references therein. Many of the models 
with preferred frame effects are similar to the types of theories that give rise to dynamical Lorentz 
violation. For example, the vector-tensor theories of Will, Hellings and Nordtvedt |278l 12371 1145) 
have many similarities to the aether theory of section l4~4l The aether model's wave spectrum has 
been calculated in [1571 1136) and limits from the absence of Cerenkov emission of these modes by 
cosmic rays has been studied in [105) (see section If). 5. 7[) . Other consequences of dynamical Lorentz 
violation in Riemann-Cartan spacetimes have been examined in [5 7) . 

Unfortunately, few constraints currently exist on dynamical Lorentz violation from gravitational 
wave observations as the spectrum is only part of the story. Currently, the expected rate of produc- 
tion of these modes from astrophysical sources as a function of the coefficients in the Lagrangian 
is unknown. However, both the energy loss from inspiral systems due to gravitational radiation 
and gravitational wave observatories such as LIGO and LISA should produce strict bounds on the 
possibility of dynamical Lorentz violating fields. 32 We note that aether type theories seem to be 
free of certain obvious problems such as a van Dam-Veltman-Zakharov type discontinuity [136] . 
The theories can therefore be made arbitrarily close to GR by tuning the coefficients to be near 
zero. 

7.2 Cosmology 

Cosmology also provides a way to test Lorentz violation. The most obvious connection is via infla- 
tion. If the number of e-foldings of inflation is high enough, then the density fluctuations respon- 
sible for the observed cosmic microwave background (CMB) spectrum have a size shorter than the 
Planck scale before inflation. It might therefore be possible for trans-Planckian physics/quantum 
gravity to influence the currently observed CMB spectrum. If Lorentz violation is present at or 
near the Planck scale (as is implicit in models that use a modified dispersion relation at high 

31 There are exceptions, for example see )151| . Here gravity is modified by a Chern-Simons form, yet there are still 
only two gravitational wave polarizations. The only modification is that the intensity of the polarizations differs 
from what would be expected in general relativity. 

32 It has also been proposed that laser interferometry may eventually be capable of direct tests of Planck suppressed 
LV dispersion 1261 . 
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energies |5T3] ). then the microwave background may still carry an imprint. 33 A number of authors 
have addressed the possible signatures of trans-Planckian physics in the CMB, for a sampling 
see [Ml EH1 EH ESI IMI HSU ISS1 E2S1 and references therein. While the possibility of such 
constraints is obviously appealing, the CMB imprint (if any) of trans-Planckian physics, much 
less Lorentz violation, is model dependent and currently the subject of much debate. 34 In short, 
although such cosmological explorations are interesting and may provide an eventual method for 
ultra-high energy tests of Lorentz invariance, for the purposes of this review we forego any more 
discussion on this approach. 

A simple low energy method to limit the coefficients in the aether model (|49|l that is less 
fraught with ambiguities has been explored by Carroll and Lim |76) . They consider a simplified 
version of the model l]49[l without the C4 term and choose the potential V{u a u a ) to be of the form 
\{u a u a — a 1 ) where A is a Lagrange multiplier. Without loss of generality, we can rescale the 
coefficients Cj in (|49|) to set a 2 = 1. In the Newtonian limit Carroll and Lim find that Newton's 
constant as measured on earth is rescaled to be 

9C 

G°n = ~ • (H5) 

In comparison, the effective cosmological Newton's constant is calculated to be 

or 1 

riobs (~\~\(\\ 

^ cosmo o ( 1 o 1 \ ' V / 

2 - (ci + 3c 2 + c 3 J 

The difference between the cosmological and Newtonian regimes implies that we have to adjust 
our measured Newton's constant before we insert it into cosmological evolution equations. Such 
an adjustment modifies the rate of expansion. A change in the expansion rate modifies big bang 
nucleosynthesis and changes the ratio of the primordial abundance of 4 He to H . By comparing 
this effect with observed nucleosynthesis limits, Carroll and Lim are able to constrain the size 
of 01,02,03. In addition to the nucleosynthesis constraint, the authors impose restrictions on the 
choice of coefficients such that in the preferred frame characterized by u a the perturbations 5u a 
have a positive definite Hamiltonian, are non-tachyonic, and propagate subluminally. With these 
assumptions Carroll and Lim find the following constraint 

< (14ci + 21c 2 + 7c 3 + 7c 4 ) < 2 (117) 

where the C4 dependence has been included for completeness. 



7.3 PPN Parameters 

Preferred frame effects, as might be expected from Lorentz violating theories, are nicely summa- 
rized in the parameterized post-Newtonian formalism, otherwise known as PPN (for a description, 
see |277) or |27fjp . The simplest setting in which the PPN parameters might be different than GR 
is in the static, spherically symmetric case. For static, spherically symmetric solutions in vector- 
tensor models the only PPN parameters that do not vanish are the Eddington-Robertson-Schiff 
(ERS) parameters 7 and f3. For GR, /3 = 7 = 1. The ERS parameters for the general Hellings- 
Nordvedt vector-tensor theory [145] are not necessarily unity |27fi) . so one might expect that the 

33 The B-mode polarization of the CMB might also carry an imprint of Lorentz violation due to modifications in 
the gravitational sector 12061 . 

34 The above approach presumes inflation and speculates about the low energy signature of Lorentz violating 
physics. Lorentz violation can also be a component in the so-called variable speed of light (VSL) cosmologies (for 
a review see |2()t) ) which are a possible alternative to inflation. Some bounds on VSL theories are known from 
Lorentz symmetry tests, but in these cases the VSL model can be equivalently expressed in one of the frameworks 
of this review. 
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constrained aether model also has non-trivial ERS parameters. However, it turns out that the 
constrained aether model with the Lagrange multiplier potential also has (3 = 7 = 1 for generic 
choices of the coefficients 103 . Therefore, at this point there is no method by which the ERS 
parameters can be used to constrain Lorentz violating theories. The ERS parameters for more 
complicated theories with higher rank Lorentz violating tensors are largely unknown. 

If we move away from spherical symmetry then more PPN parameters become important, in 
particular ot\, a.%, a 3 which give preferred frame effects. In |I33| a 2 has been calculated to be 

Q'2 = -=-7 r 1 7 [2c 1 +4c2( C2 + C3 ) + c 2(3c 2 + 5c 3 + 3c4) (118) 

ZCi(Ci + C 2 + C 3 J 

+ Cl ((6c 3 - c 4 )(c 3 + c 4 ) + c 2 (6c 3 + c 4 ))] . (119) 

The observational limit on a 2 is |a 2 | < 4 • I0 -7 |277| . Barring cancellations this translates to a 
very strong bound of order I0 -7 on the coefficients Ci in the aether action. 
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8 Conclusions and prospects 



As we have seen, over the last decade or two a tremendous amount of progress has been made in 
tests of Lorentz invariance. Currently, we have no experimental evidence that Lorentz symmetry 
is not an exact symmetry in nature. The only not fully understood experiments where Lorentz 
violation might play a role is in the (possible) absence of the GZK cutoff and the LSND anomaly. 
New experiments such as AUGER, a cosmic ray telescope, and MiniBooNE |112) . a neutrino 
oscillation experiment specifically designed to test the LSND result, may resolve the experimental 
status of both systems and allow us to determine if Lorentz violation plays a role. 

Terrestrial experiments will continue to improve. Cold anti-hydrogen can now be produced 
in enough quantities [1181 12 7) for hydrogen/anti-hydrogen spectroscopy to be performed. The 
frequency of various atomic transitions (1S-2S, 2S-nd, etc.) can be observationally determined 
with enough precision to improve bounds on various mSME parameters [fill 1256) . Spectroscopy 
of hydrogen-deuterium molecules might lead to limits on electron mSME parameters an order of 
magnitude better than current cavity experiments |228| . 

There are proposals for space based experiments (c.f. |59l 1195) ) that will extend current 
constraints from terrestrial experiments. Space based experiments are ideal for testing Lorentz 
violation. They can be better isolated from contaminating effects like seismic noise. In a micro- 
gravity environment interferometers can run for much longer periods of time as the cooled atoms 
in the system will not fall out of the interferometer. As well, the rate of rotation can be controlled. 
Sidereal variation experiments look for time dependent effects due to rotations. In space, the rate 
of rotation can be better controlled, which allows the frequency of any possible time-dependent 
signal to be tuned to achieve the best signal-to-noise ratio. Furthermore, space based experiments 
allow for cavity and atomic clock comparison measurements to be combined with time dilation 
experiments (as proposed in OPTIS |195) ). thereby testing all the fundamental assumptions of 
special relativity. The estimated level of improvement from a space based mission such as OPTIS 
over the corresponding terrestrial experiments is a few orders of magnitude. 

Another possibility for seeing a novel signal of Lorentz violation is in GLAST 260 . GLAST 
is a gamma ray telescope that is very sensitive to extremely high energy gamma ray bursts. As 
we have mentioned, DSR evades almost all known high energy tests of Lorentz invariance. If the 
theoretical issues are straightened out and DSR does eventually predict a time of flight effect then 
GLAST may be able to see it for some burst events. An unambiguous frequency to time-of-arrival 
correlation linearly suppressed in the Planck energy, coupled with the observed lack of birefringence 
at the same order, will be a smoking gun for DSR, as other constraints forbid such a construction 
in effective field theory \2'M)\ . 

The question that must be asked at this juncture in regards to Lorentz invariance is: when 
have we tested enough? We currently have bounds on Lorentz violation strong enough that there 
is no easy way to put Lorentz violating operators of dimension < 6 coming solely from Planck scale 
physics into our field theories. It therefore seems hard to believe that Lorentz invariance could be 
violated in a simple way. If we are fortunate, the strong constraints we currently have will force us 
to restrict the classes of quantum gravity theories/spacetime models we should consider. Without 
a positive signal of Lorentz violation, this is all that can reasonably be hoped for. 
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